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Preface 



This volume has grown from a conference entitled Harmonic Maps, Minimal Sur- 
faces and Geometric Flows which was held at the Universite de Bretagne Occi- 
dentale from July 7th-12th, 2002, in the town of Brest in Brittany, France. We 
welcomed many distinguished mathematicians from around the world and a dy- 
namic meeting took place, with many fruitful exchanges of ideas. 

In order to produce a work that would have lasting value to the mathematical 
community, the organisers decided to invite a small number of participants to 
write in-depth articles around a common theme. These articles provide a balance 
between introductory surveys and ones that present the newest results that lie at 
the frontiers of research. We thank these mathematicians, all experts in their field, 
for their contributions. 

Such meetings depend on the support of national organisations and the local 
community and we would like to thank the following: the Ministere de VEducation 
Nationale, Ministere des Affaires Etrangeres, Centre National de Recherche Scien- 
tifique (CNRS), Conseil Regional de Bretagne, Conseil General du Finistere, Com- 
munaute Urbaine de Brest, Universite de Bretagne Occidentale (UBO), Faculte des 
Sciences de VUBO, Laboratoire de Mathematiques de VUBO and the Departement 
de Mathematiques de VUBO . Their support was generous and ensured the success 
of the meeting. 

We would also like to thank the members of the scientific committee for 
their advice and for their participation in the conception and composition of this 
volume: Pierre Berard, Jean-Pierre Bourguignon, Frederic Helein, Seiki Nishikawa 
and Franz Pedit. 

A large conference requires an enormous amount of work behind the scenes; 
we thank Annick Nicolle , secretary for the Laboratoire de Mathematiques at the 
UBO, for her administrative skills which ensured the smooth running of the meet- 
ing. 

Finally, more than ninety mathematicians participated at the conference and 
we hope also that this volume will be a tribute to their enthusiasm. 

The organisers Paul Baird 

Ahmad El Soufi 
Ali Fardoun 
Rachid Regbaoui 
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Introduction: Best Maps and Best Metrics 

Paul Baird, Ahmad El Soufi, Ali Fardoun and Rachid Regbaoui 



1. Introduction 

Among the most beautiful problems of mathematics are those which we can visu- 
alise, where the solution contains the faith of our intuition, and yet where achieving 
this solution proves tantalisingly elusive. Such problems demand the introduction 
of new techniques and ideas in order to advance. The uniformisation of geomet- 
ric structures has provided such an impetus for the development of differential 
geometry over many years. We refer back to Riemann and the problem of recog- 
nising Riemann surfaces and the conformal structures defined upon them; to the 
influence of physics leading to the notion of Einstein metric ; to the problem of 
identifying geometric structures on 3-manifolds and that of finding a conformal 
metric of constant scalar curvature - the so-called Yamabe problem. Naturally, in 
looking for “best” geometric structures, one is lead to look for “best” maps as well, 
for example harmonic maps and p-harmonic maps , eigenvalue problems or more 
general solutions of elliptic equations. The unifying feature of all these problems 
is their variational characterization. Solutions arise as extremals of an “action” 
integral, perhaps with constraints. 

The articles in this volume are devoted to the study of such problems; some 
papers provide an overview whereas others present new work and lead the reader 
to the very latest developments. 



2. Best metrics 

The basic object is a manifold M m ; unless we specify otherwise, this will be a 
smooth connected Hausdorff paracompact manifold without boundary. By taking a 
partition of unity, we can endow M m with a smooth metric g. We will be concerned 
with metrics that are positive definite , i.e., g(v,v) > 0 for every non-zero tangent 
vector v. Our interest is to find relations between the topology of M m and its 
geometry. Here intuition plays a role: we can “see” that a sphere is fundamentally 
different to a torus, both topologically (there are curves on the torus that cannot 
be shrunk to a point) and geometrically (the standard sphere “bends” in the same 
way everywhere, whereas the torus does not). To make these intuitions precise, we 
introduce the notion of curvature. 
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If V denotes the Levi-Civita connection on (M m ,g), then the Riemannian 
curvature is the 3-covariant 1-contravariant tensor field R M given by 

R m (X, Y)Z = VxVyZ - VyVx^ - V[X,Y]Z (X, Y,Z e C°°TM). 

If X, Y are orthonormal, the sectional curvature of the plane spanned by X and 
Y is the number g(R M (X, Y)Y, X). 

The Ricci tensor is the symmetric tensor field Ric M defined by 

rri 

Ric M (X,Y) = Tv g g(R M (X,-)-,Y) = Y) , 

2=1 

where {e*} is an orthonormal frame. Finally, the scalar curvature is the function 
S M = Scal M defined by 

m 

Scal M = Tr s Ric M = ^Ric^e*) . 

2=1 

For a surface, we have the basic identity: 

R M (X,Y)Z = K M (g(Y,Z)X-g(X,Z)Y) , 

where K M denotes the Gaussian curvature of M 2 (the same formula holds in ar- 
bitrary dimension, when the manifold M has constant sectional curvature K M ). 
Taking traces, Ric M = K M g and Scal M = 2K M , so there is essentially just one 
curvature for a surface, the Gaussian curvature. (Here we are restricting our discus- 
sion to instrinsic curvature; of course for an immersed manifold we could discuss 
the mean curvature , for example.) 

In dimension 2, we have a clear idea of what we mean by “best” metric - one 
of constant Gaussian curvature. The Riemann mapping theorem (or uniformisa- 
tion theorem) asserts that every simply-connected domain of the complex number 
sphere S 2 with a minimum of two boundary points can be mapped conformally 
and diffeomorphically onto the interior of the unit disc (see [21] for a proof of this 
- apparently Riemann gave no rigorous proof!). The consequences of this theorem 
are far reaching. For example, every compact Riemann surface M 2 , is conformally 
equivalent to E/F, where E = S' 2 , R 2 ,# 2 is one of the two-dimensional space 
forms and T is a discrete group of isometries acting properly and discontinuously. 
This gives us a complete picture of a compact Riemann surface and implies that 
each one can be deformed in its conformal class to one of constant curvature. 

The problem of deforming a metric on a surface to one of constant Gaussian 
curvature has a variational formulation. Indeed, if (M 2 ,go) is a compact Riemann 
surface without boundary with scalar curvature So, then a conformal deformation 
of the form g — e 2u go produces a metric with scalar curvature S = e~ 2u (— 2Aou + 
So), where A 0 is the Laplace-Beltrami operator with respect to the metric go- 
Thus, in order to find a metric of constant curvature, we are required to solve the 
partial differential equation 



— 2Aq u + Sq = e 2w . 
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The left-hand side is now the gradient of the functional 

I(u) = / \Vu\ld/i 0 + 2 S 0 udg 0 , 

JM JM 

subject to the constraint f M e 2u dg 0 = So f M dgo, where dgo denotes the canonical 
measure and Vu the gradient of u with respect to the metric go (see [22], [2]). 
Hamilton gave an elegant way of obtaining a solution from the evolution equation 

! = (r-S)9, (1) 

where r is the average value of S [10]; we shall discuss such equations in more 
detail below. 

In higher dimensions, it is therefore natural to formulate the following prob- 
lem: given a compact Riemannian manifold (M m ,go), find a metric g — u 4 S n ~ 2 )g 0 
with constant scalar curvature. This problem was first posed by Yamabe in 1960 
[26] and is now known as the Yamabe problem. Yamabe’s work had an error noticed 
by Trudinger [24], who solved the problem in a special case. Since that time there 
have been contributions from a number of authors. It was solved completely in the 
compact case by Aubin [1] and Schoen [19]. The delicate aspect of this problem 
concerns the critical exponent in the Sobolev embedding theorem, as follows. 

The equation to solve is now 

-A 0 U + -f—YSoU = A u (n+2)/(n- 2 ) ( 2) 

4 (n — 1) 

giving a metric of constant scalar curvature S = 4 (n — l)A/(n — 2). This has the 
following variational characterisation. 

Let H 2 (M) be the Sobolev space of functions in L 2 whose first derivatives 
are also in L 2 . Let I be the functional associated to equation (2): 

I(u) = [ |Vu|od/i 0 + -Jf — ~tt [ Sou 2 dgo • 

Jm 4(n — 1) J M 

By the Rellich-Kondrakov theorem, the embedding H 2 (M) ^ L q (M ), q G (2, ^ 2 ) 
is compact; which leads to a solution of the subcritical equation 

nn O 

-A ou + — 7 tSou = X q u q ~ l , 

4 (n — 1) 

subject to the constraint u G H q = {u G H li2 (M) : f M u q dfio = 1} , where X q = 
inf {I(u) : u e H q }. Subtle arguments are required to extend this to the critical 
exponent q = 2n/(n — 2). Essentially, one is required to show that there is a smooth 
everywhere positive function u for which the Yamabe invariant go = inf U eu H u ) 
is achieved, where H = {u G Hf(M) : f M \u\ 2n ^ n ~ 2 ^dgo = 1} - see [12] for an 
excellent account. Recent developments concerning critical Sobolev exponents are 
discussed in the article of E. Hebey in this volume [13]. 

Although a beautiful problem in dimension m > 2, the solution to the Yam- 
abe problem does not give us a conclusive description of manifolds. For example 
the manifolds S 3 and S 2 xR both have constant positive scalar curvature. By fac- 
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toring through discrete groups of isometries we can obtain compact 3-manifolds 
which are fundamentally different. From the first, we could obtain the Lens spaces, 
for example, which cannot be obtained from the second (see [20] for an account of 
such matters). 

In order to look for “best” metrics in higher dimensions, it is natural to 
consider the functional: 

K{g) = f Scal M dn(g) (3) 

J M rn 

defined on a suitable space of metrics. In order to avoid trivial solutions, we need 
to impose the constraint f Mm d ju(g) = 1 , i.e., that the volume remain constant. 
On taking a variation of (3) with respect to the metric g subject to the constraint, 
we obtain the Euler-Lagrange equations in the following form: 

^Scal M g — Ric M = ag 

where a is a constant. Solutions therefore occur when Ric(g) = A g, for some 
function A. In dimension m > 2, the Bianchi identities now imply that A must be 
constant and we are therefore led to the equation 

Ric {g) = eg 

with c constant. A metric satisfying this equation is called an Einstein metric. Un- 
fortunately in dimension m = 3, this condition does not give us enough flexibility 
to pick out all compact manifolds, since it implies that M 3 must have constant 
sectional curvature which in turn implies that M 3 = E/T, where E = S 3 ,R 3 ,iL 3 
is one of the three-dimensional space-forms and T is a discrete group of isome- 
tries acting discretely and properly discontinuously. However, in dimension 4, the 
construction of Einstein metrics is a rich field of study, with much recent progress. 

In dimension 3 the picture is rather delicate and we are guided by the con- 
jecture of Thurston. A model geometry is a smooth simply connected manifold 
E together with a Lie group G of diffeomorphisms of E which act transitively 
on E with compact point stabilisers, such that G is maximal in the sense that 
it is not contained in any larger group of diffeomorphisms of E with compact 
point stabilisers. It is also a requirement that there exists at least one compact 
quotient, i.e., there exists a subgroup H of G such that E/H is a compact man- 
ifold. Two such geometries are considered equivalent if there exists a diffeomor- 
phism between them which intertwines the group actions. Thurston has classified 
the geometries; up to equivalence there are eight of them and they are given 

by S' 3 , R 3 , iL 3 , S 2 x R, H 2 x R, Nil, SL 2 (R), Sol. The geometrization conjecture of 
Thurston asserts that any compact 3-manifold is made up of pieces, separated by 
2-spheres or 2-tori, each piece covered by one of the eight geometries. The resolu- 
tion of this conjecture in the affirmative would imply the Poincare conjecture. 

A powerful result, in the spirit of the above ideas, was proved by R. Hamilton 
[9]. He established that if (M 3 ,g) is compact with Ric(g) > 0, then the metric g 
can be smoothly deformed to one of constant sectional curvature. This means that 
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M 3 is diffeomorphic to 5 3 /T, where T is a discrete group of isometries acting 
properly discontinuously. His method was to allow the metric to evolve using the 
so-called Ricci flow: 

% = -2Ric(flr) + -rg (n = 3) , 
ot n 

where r = f M S d/i / f M d/i is the average scalar curvature (this in fact restricts to 
the evolution equation (1), in the case of dimension 2). Since the right-hand side 
is a semi-elliptic operator in the components of the metric, this equation should 
be viewed as a kind of heat equation . Intuition suggests that in the limit, we will 
obtain a metric of constant Ricci curvature. As already noted, in dimension 3, this 
means that the manifold has constant sectional curvature and the result follows. 
More recently, in order to attack Thurston’s geometrisation conjecture, Hamilton 
has proposed a program involving the unormalised Ricci flow: 

i=- 2Ri ^»- 

The article of D. Knopf [15] in this volume, discusses the very latest work on this 
method and notably the possible resolution of the geometrisation conjecture by 
G. Perelman. 

Geometrically, fixed points (up to diffeomorphism and dilation) of the Ricci 
flow are interesting objects (see [15]); they are metrics g which satisfy the equation 

— 2Ric(g) = C x g + ag , 

for a constant a and a vector field X. The article of Gastel [8] shows how to obtain 
examples of solitons by a warped product construction. 



3. Best maps 

In seeking “best” maps, we should note that they may be defined in terms of a 
metric and therefore a “best” metric may play a role - i.e., the two problems are 
interrelated, and indeed we find that their variational characterisation and descrip- 
tion as solutions to an elliptic partial differential equation have many analogies. 

Let us briefly recall the main theorem of Hodge theory (see, for example 
[25]). Let (M m ,g) be a compact Riemannian manifold; then each element of the 
kth order de Rham cohomology group H k (M, R) contains a unique harmonic 
representative. For a differential form uj G H k (M, R), the Laplacian is defined 
by A uj = dd*cj + d*dcj, where d* is the codifferential operator characterized by 
f M ( d*cu , rj) d ix(g) = f M ( uj , d rj) d g(g) for all (k - l)-forms rj. A harmonic form to is 
one which satisfies the elliptic equation A uj = 0, equivalently, using the compact- 
ness of M, harmonic forms are characterised by the equations do; = d*a; = 0. This 
remarquable theorem gives a direct link between the topology of M and solutions 
to an elliptic differential equation. 

Harmonic mappings were introduced in 1954 by Sampson in an unpublished 
MIT report, in an attempt to obtain a homotopy version of Hodge theory. Another 
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definition was proposed by J. Nash and H.E. Rauch (unpublished). F.B. Fuller 
(1954) gave the same definition as Sampson and obtained the first variation for- 
mula [7]. 

For a map p : (M m ,g) — » (A^ n , h) between Riemannian manifolds, we define 
the energy density by e(p) = ||d(/?| 2 , where \dp\ 2 = Tr g p*h is the Hilbert-Schmidt 
norm of dip. If M is compact, we can now define the (total) energy to be the integral 

E{v) = [ e(<p)dfj,(g) . 

J M 

The Euler-Lagrange equations associated to this functional are 

— d*d(/? := r(p) = 0 . 

We note that, since automatically dd(^ = 0, on a compact manifold, the Euler- 
Lagrange equations are equivalent to Ad p = 0, where A is now the Laplacian 
acting on vector bundle valued forms (viewing dip £ C°°(T*M <S> p~ l TN)). By 
analogy with Hodge theory, we can formulate the fundamental problem of har- 
monic maps as follows. Given a homotopy class Ti £ C°°(M, N) of smooth maps 
from (M,g) to (N, h), find a harmonic representative p £ hi (We avoid discussion 
of the finer points of regularity theory in this introduction and suppose all maps 
smooth.) 

The first major study of harmonic maps between Riemannian manifolds was 
made by Eells and Sampson in 1964 [5]. Their most striking theorem was the 
following: Let (M m ,g) and ( N n , h) be compact Riemannian manifolds and suppose 
that (N, h) has non-positive sectional curvature. Then, for any homotopy class 
H £ C°°(M, N), there is a harmonic representative p ^TL. 

As important as the result, was the method of proof. This was to use the 
heat equation: 

Tt= T{lf) - (4) 

The technique encapsulates the intuition referred to in the opening paragraph. We 
may imagine a closed curve sliding down the tube of a trombone until, at time 
infinity, it arrives at an absolute minimum of the energy where the tube is most 
narrow. The resulting map will be a geodesic - a special case of a harmonic map. 
In a similar way in which solitons occur (self-similar solutions), “bubbles” can 
occur in the harmonic map flow. This is studied in the article of Topping [23] and 
appears in a more general context in the contribution of Hebey [13]. The natural 
objects with which to describe this type of phenomena are called rectifiable scans. 
The setting is now the domain of geometric measure theory, and the application 
of these to Plateau-type minimisation problems is discussed in the contribution of 
De Pauw and Hardt [3]. 

It is known that some homotopy classes of maps cannot contain a harmonic 
representative, for example, there can be no harmonic map of degree 1 from the 2- 
torus to the 2-sphere [6]. However, a striking result of Eells and Ferreira, established 
in 1991 [4] realises Sampson’s original goal, provided we allow deformations of the 
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metric on the domain (so reinforcing our idea of interdependence between best 
metrics and best maps): given two compact Riemannian manifolds ( M,g ) and 
(TV, h ) with dim M > 3 and a homotopy class hi of smooth maps between them, 
there is a smooth metric g conformally equivalent to g and a map (p G H such 
that cp : (M, g) — » (TV, h) is harmonic. Their method of proof was to consider the 
modified energy functional 

\ f (i + \M 2 ) p/2d Kg) ■ (5) 

z J M 

Now, if p > 2 dim M, general principles from the direct method in the calculus of 
variations show that this attains its minimum in each homotopy class. 

The appearance of the functional (5) illustrates the importance of the study 
of more general types of functionals. The most natural generalisation of the energy 
of a map, is the p- energy, defined by 

E P {<f) = - f \dip\ p dn(g) . 

P JM 

We may for example hope that special properties of harmonic maps in dimension 2 
may be transferred to dimension p for critical points of this functional (p-harmonic 
maps). However, the Euler-Lagrange equations r p ((p) = d*(|d(/?| p-2 d(/?) = 0 degen- 
erate when d(p = 0 and the existence and regularity properties of p-harmonic maps 
becomes a difficult domain of study. The article of Hungerbiihler [14] in this vol- 
ume investigates the existence and regularity of p-harmonic maps by means of the 
heat flow. 

Other aspects of harmonic maps, notably harmonic maps in complex Finsler 
geometry and harmonic maps from a flat complex are discussed in the articles of 
S. Nishikawa [18] and C. Mese [17], respectively, in this volume. 

Finally, in the spirit of the evolution methods discussed above, one of the most 
successful programs in recent years has been the evolution of a hypersurface by its 
mean curvature, the so-called mean curvature flow. The article by C. Mantegazza 
[16] discusses the evolution of networks in the plane. This is an extension of these 
ideas in a new and imaginative way and represents the most recent work in this 
area. 

The papers in this volume indicate the progress that has been made in recent 
years towards the goal of understanding the interplay between the topology and 
the geometry of manifolds. They are written by experts in their respective fields 
and we hope will provide inspiration and a platform for others to tackle some of 
the many beautiful problems that still remain. 
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Bubbles over Bubbles: A C°-theory 
for the Blow-up of Second Order Elliptic 
Equations of Critical Sobolev Growth 

Emmanuel Hebey 



Abstract. Let (M, g ) be a smooth compact Riemannian manifold of dimension 
n > 3, and A^ = — div^V be the Laplace-Beltrami operator. Let also 2* be 
the critical Sobolev exponent for the embedding of the Sobolev space Hf (M) 
into Lebesgue’s spaces, and h be a smooth function on M. Elliptic equations 
of critical Sobolev growth like 

A g u + hu = u 2 -1 

have been the target of investigation for decades. A very nice ifi-theory for 
the asymptotic behaviour of solutions of such an equation is available since the 
1980’s. We discuss here the C°-theory recently developed by Druet, Hebey and 
Robert. Such a theory provides sharp pointwise estimates for the asymptotic 
behaviour of solutions of the above equation. 



Let (M, g ) be a smooth compact Riemannian manifold of dimension n > 3. 
We denote by H 2 (M) the standard Sobolev space of order two for integration and 
order one for differentiation. We let (h a ) be a sequence of C 0,e functions on M, 
0 < 6 < 1, and consider equations like 

A g u + h a u = u 2 -1 (1) 

where A^ = — div^V is the Laplace-Beltrami operator, 2* = 2 n/(n — 2) is the 
critical Sobolev exponent for the embedding of the Sobolev space H 2 (M) into 
Lebesgue’s spaces, and u is required to be positive. We regard (1) as a possible 
model for more general second order elliptic equations of critical Sobolev growth. 
We let (u a ) be a bounded sequence in Hf(M) of solutions of (1) in the sense that 
for any a, 

A g u a + h a u a = u 2 a -1 

and ||ua||//2 < A where A > 0 is independent of a. We also assume that the h^s 
are uniformly bounded and that they converge in L 2 to some limiting function hoc. 
Then, thanks to Struwe [32], we know how to describe the asymptotic behaviour 
of the s in H 2 (M) as a — > +oo. More precisely, it follows from Struwe [32] that, 
up to a subsequence, 



u a =U 0 + J2 B a + R t 



( 2 ) 
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where u° is a solution of the limit equation 

A g u + hooU = u 2 ~ 1 

the sum in the right-hand side of (2) is a finite sum over i, B l a is a bubble, 
obtained by rescaling fundamental positive solutions of the Euclidean equation 
A u = u 2 _1 , and the R a s are lower order terms which converge strongly to 0 in 
Hf(M). This asymptotic description provides a very satisfactory JT^-theory for the 
asymptotic behaviour of solutions of equations like (1). Let us assume now that the 
s converge C 0 ' 6 to hoo for some 0 < 6 < 1. An important issue in the study of 
equations like (1) is to get a theory in which the above asymptotic description holds 
also in the C°-space, where pointwise estimates are involved. Such a C°-theory was 
recently developed in Druet, Hebey and Robert [12, 13]. We discuss this theory 
in these notes. Important applications of the theory are when dealing with sharp 
Sobolev inequalities. Monographs on sharp constant problems are Druet-Hebey 
[10], and Hebey [17]. Other possible applications are when dealing with the energy 
function, see for instance Hebey [18], or when dealing with compactness results as 
in Schoen [28, 29]. Very nice applications of the theory in such directions are in 
Druet [9]. 



1. Background material for the Euclidean space 

We assume that n > 3. Thanks to the Sobolev embedding theorem [31], we refer 
also to Gagliardo [15] and Nirenberg [23], there exists K > 0 such that for any 
smooth function u with compact support in R n , 

N|2* <K\\Vu\\ 2 (3) 

where 2* = 2 n/(n — 2) is as above. Another classical result is that if K n stands 
for the sharp constant in (3), then 



Y n(n - 2 )uj 2 r! n 

where c o n is the volume of the unit n-sphere. Possible references on sharp Sobolev 
inequalities, including the Riemannian case, are Druet-Hebey [10] and Hebey [17]. 
The sharp Euclidean Sobolev inequality reads as 

\\u\\ 2 * <K n \\Vu\\ 2 . 



The extremal functions for this inequality are known. Let u be the function given 
by 



n — 2 



2 



1 + 



l*l a 

n(n — 2) 



u(x) = 



(4) 
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Then, see for instance Bliss [2], u is an extremal function for the sharp Euclidean 
Sobolev inequality. Moreover, as it was shown by Caffarelli, Gidas and Spruck [6], 
see also Obata [24] , u is the only positive solution of the critical Euclidean equation 

A u — u 2 -1 (5) 

which is such that u(0) = max^n u = 1. All the positive solutions of (5) are then 
given by 

u a( x ) = u(\(x — a)) (6) 

where A is any positive real number, a is any point in the Euclidean space, and u is 
the fundamental solution of (5) given by (4). In other words, positive solutions of 
(5) come from rescalings of the fundamental solution u given by (4). As a remark, 
the energy E(u) — ||u|| 2 * of u is given by E(u) = 



2. Background material for the unit sphere 



The case of the unit sphere is closely related to the case of the Euclidean space. 
Let (*S n , h) be the unit n-sphere, n > 3. Let Sh be the scalar curvature of h. Then 
Sh is constant and Sh = n(n — 1). We consider the conformal scalar curvature 
equation 



A hU + 



n(n — 2) 9 * 

: ~ u — u 

4 



-l 



so that, with respect to (1), 



k _ n-2 0 

for all a. A classical result is that positive solutions of this equation are known. 
Let xo G S n and r be the distance to xq. Then, for any (3 > 1, 

Up = ( n ^ n 4 2 \ p 2 - !)") (/3-cosr) 1 ~ 1 



is a solution of the above equation. The energy E(u) = \\u\\ 2 * of up is given by 



E(up) = Kp n - 2 ^ 2 



where K n is the sharp constant as in the preceding section. A possible refer- 
ence in book form on such topics is Hebey [17]. It is easily seen that up — > 0 in 
Cfoc (S n \{xo}) as (3 — > 1. On the other hand, up(x o) — > +oo as (3 1. Thus the 

up s develop a singularity at xo as (3 — > 1. We say that the up's blow-up at xo- Let 
xp and (ip be given by 

up(xp) = ma xup(x) = /A 2 . 
x£S n P 
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Then xp — xq and (Xp — y n ^^ 2 ){p+i) • particular, pp — > 0 as (3 — > 1. Let 5/3 be 
the function given by 



Bp(x) 



M/3 



m| + 



d h (x@,x) 2 
n(n— 2 ) 



n-2 

2 



We refer to 5/3, see the following sections, as the standard bubble with respect to 
xp and [ip. It is easily seen in this particular situation that there exists C > 1 
such that 

-^Bp(x) < up{x) < CBp(x) 

for all x G S n and all (3 > 1. Similarly, for any e > 0, there exists > 0 such that 

< u p(x) < (1 + e)Bp(x) 

for all x G B Xo (8 e ) when (3 > 1 is close to l.Let Rp be the function given by 



up = Bp T Rp . 



Then, another easy claim in this particular situation is that 

Rp^O in Hf(S n ) 

as (3 — > 1. On the other hand, in general, Rp 0 in C°(S n ) as (3 — ► 1, and the 
5 / 3 ’s are not even bounded in L°°(S n ). Let us assume for instance that n > 7, 
and let(yp) be a sequence of points in S n such that dh(xo,yp) 2 = ((3 — l) 1+£ 
where ,e > 0 is such that e < Then, Rp{yp) — > +oo as /? — > 1. A more 

complete information is as follows. When n > 7, the Rp's are not bounded in 
L°°(S n ). When n = 6, the 5/3 ’s are bounded in L°°(5 n ), but they do not converge 
to 0 in C°(S n ) as (3 — » 1. When n = 3,4,5, the Rp ’s converge to 0 in C 0 ^ 72 ) 
as (3 — ■» 1. Other examples of blowing-up sequences are in Druet-Hebey [11] and 
Druet- Hebey- Robert [12]. 



3. Blow-up theory in the Hf - Sobolev space 

When dealing with Sobolev spaces, we do not really need the notion of exact 
solutions. The more general notion of a Palais- Smale sequence works also. Let 
(M,g) be a smooth compact Riemannian manifold of dimension n > 3, and (h a ) 
be a sequence of continuous functions on M. We consider equation (1), and assume 
that there exists C > 0, and a continuous function hoo on M such that the following 
holds: 

(Al) for any a, and any x G M, \h a (x)\ < C, and 
(A2) h a — > hoo in L 2 (M) 

as a — > Too. Then, of course, it follows from (Al) and (A2) that the h a ’s converge 
to in all the Lebesgue spaces L P (M), p > 1. In other words, we assume that the 
h a ^ s are bounded in L°°(M) and that they converge in Lebesgue spaces to some 
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limiting function . Let u be a function in Hi(M). We define the functional Ig 
by 

I g( u ) = l[ (|Vm| 2 + h a u 2 ) dv g - L f \ufdVg 
z J M Z J M 

so that Ig is the formal primitive of equation (1). We let ( u a ) be a sequence 
of functions in the Sobolev space By definition, we say that ( u a ) is a 

Palais-Smale sequence for equation (1) if the two following conditions hold: 

(i) Ig(u a ) is bounded with respect to a, and 

(ii) DIg(u a ) -> 0 strongly 

as a — ► Too. This well-known notion of a Palais-Smale sequence extends the notion 
of a sequence of solutions of equation (1), since for an exact solution, the differential 
of Ig is zero. The goal in this section is to characterize the behaviour as a — > +oo 
of such Palais-Smale sequences. In order to do that, we need the important notion 
of a bubble. We let r]s be a cut-off function in M n such that rjs{x) = 1 if \x\ < 5 and 
r]s(x) = 0 if |x| > 25, where 0 < 5 < i g /2 and i g is the injectivity radius of (M,g). 
We let (x a ) be a sequence of points in M, and (R a ) be a sequence of positive real 
numbers such that +oo as a -> +oo. We let also u % be a positive solution 

of the critical Euclidean equation (5), so that u * is given by (6). Then we define 
a bubble as a sequence (B a ) of functions given by an equation like 

n — 2 

B a (x) = R^~ rj 6iXa (x)u^(R a exp~^(x)) 

where exp Xa is the exponential map at x a , and r]$ :Xa (x) = rjs(exp~^(x)) . In par- 
ticular, since we do know what u £ is, we do know the exact shape of bubbles. 
Returning to our question, we want to characterize the behaviour as a — » +oo of 
a Palais-Smale sequence for equation (1). We let 

A g u + hooU = u 2 -1 (7) 

be the limit equation we get when we replace h a by h ^ in (1). The result which 
answers our question is as follows. Related references are Brezis-Coron [3, 4], Lions 
[22], Sacks-Uhlenbeck [27], Schoen [29], and Wente [33]. 

Theorem 3.1 (Struwe, [32]). Assume (Al) and (A2). Let (u a ) be a Palais-Smale 
sequence for ( E a ), with u a > 0 for all a. Then there exists some integer k, there 
exists a nonnegative solution u° of the limit equation (7), and there exist k bubbles 
(B l a ), i running from 1 to k, such that , up to a subsequence, 

k 

u a = u° + B l a + R a 
2=1 

for all a, where R a — > 0 in (M) as a +cxd. 

According to this theorem, a Palais-Smale sequence for (1) splits into a so- 
lution of the limit equation, which may be zero, a sum of bubbles that come 
from solutions of the critical. Euclidean equation (5), and lower order terms in the 
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Hy -Sobolev space. As an important remark, the energy splits also. Another formu- 
lation of this is that bubbles do not interact at the H\ -level. Then, noting that the 
energy of a bubble is somehow a constant independent of the bubble, we do get a 
classification of the energy levels for which compactness does not hold. This gives 
in turn conditions under which u° is not the zero function, and thus, conditions 
under which the limit equation (7) possesses nontrivial solutions. The key point in 
the proof of this result is that given a Palais-Smale sequence, there exists a bubble 
such that the following holds. Namely that the Palais-Smale sequence minus the 
bubble is a Palais-Smale sequence whose energy is, up to lower order terms, the en- 
ergy of the original Palais-Smale sequence minus some fixed-dimensional constant 
corresponding to the energy of a bubble. Then, to some extent, k in the theorem is 
the number of bubbles we have to substract to the original Palais-Smale sequence 
in order to get a Palais-Smale sequence of small energy. Noting that we do have 
compactness when the energy is sufficiently small, this proves the theorem. 

Theorem 3.1 provides a very nice and complete H \ -theory for the blow-up 
of an equation like (1). The question now is to develop a C°-theory involving 
pointwise estimates. As already mentioned in the introduction, such a C°-theory 
was recently developed in Druet, Hebey and Robert [12, 13]. This is the subject of 
the following section. The naive idea here is that exact solutions of an elliptic PDE 
should satisfy C°-estimates and not only Hf-e stimates. Among others, a difficulty 
we have to face is that bubbles will interact at the C°-level. 



4. Blow-up theory in the C°-space - the result 

We assume in what follows that there exists h 0 0 E such that the operator 

A^ T hoo is coercive, and such that 

(A3) h a -► fioo in C 0,e (M) 

as a — > -hoo. By definition, A^ + /loo is coercive if its energy controls the H \- 
norm. This is automatically the case if is a positive function. We let (u a ) be a 
sequence of solutions of equation (1) in the sense that for any a, u a is a solution 
of (1). We assume that the u a 's are bounded in the Sobolev space Hf (M). Then, 
since Hf(M) is a reflexive space, there exists some function u° in Hf(M) such 
that, up to a subsequence, 

u a u° weakly in (8) 

as a — > -hoo. Clearly, this is easy to check, u° is a solution of the limit equation 

A g u + h^u = u r ~ x . (9) 

Another simple remark we get from standard elliptic theory is that either u a — > u° 
in C 2,d (M), in which case we do have our asymptotic description, or we do have 
that 

(A4) maxn Q — * Too 
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as a — > -boo. In this case, blow-up occurs, and this is the difficult situation in 
which we have to characterize the asymptotic behaviour of the u a ’ s. The theorem 
which answers our question, providing the C°-theory we mentioned above, is as 
follows. 



Theorem 4.1 (Druet-Hebey-Robert, [12, 13]). Let ( u a ) be a sequence of solutions 
of (1), u a > 0. We assume that the u a ’s are bounded in (M), and that (A3) 
and (A4) hold. Then there exist k G N*, converging sequences (x^ a ) in M , and 
sequences (/i 2 , a ) of positive real numbers converging to 0, i = 1, . . . , k, such that, 
up to a subsequence, 



(1 -£ a )u°(x) + 

i= 1 



Pi, a 

2 | d g (xj, a ,x) 2 

' n(n — 2) 



-2 

2 



k 

< u a (x) < (1 + e a ) u°(x) + C 

i— 1 



Pi, a 






+ 



d g {xi, a ,x ) 2 
n(n — 2) 



-2 

2 



for all x G M and all a, where u° is the solution of the limit equation (9) given by 
(8), C > 1 is independent of a and x, and (e a ), independent of x, is a sequence 
of positive real numbers converging to 0 as a — >• -boo. In particular, the u a ’s are 
C° -controlled, on both sides, by u° and standard bubbles. 



A complement to Theorem 4.1 is that C can be chosen as close as we want 
to 1 if we restrict the equation in Theorem 0.1 to small neighbourhoods of the 
geometrical blow-up points, defined as the limits of the x^ a ’s. For instance, if 
u° ^ 0, or if the uT s just have one geometrical blow-up point, then for any 5 > 0, 
there exists S £ > 0 such that, up to a subsequence, 



(1 -£ a )u°(x) + 



1 

1 + £ 



E 



Pi, a 

..2 , d g (x z ^,x) 2 

Pi, a ' n(n — 2) 



-2 

2 



k 

< U a (x) < (1 + £ a ) U°(x) + (1 + s) 

2=1 



Pi,a 

2 , d g {x i>a ,x ) 2 

Pi, a. ' n(n — 2) 



-2 

2 



for all a, all xo € 5, and all x G F? Xo (5 e ), where S is the set consisting of the 
limits of the x^ a ’s as a — > -boo, and B Xo (<S e ) is the geodesic ball of center xq and 
radius S £ . Outside the B x (5 e y s, x G <S, the ii a ’s converge C 2,0 to -u 0 . The estimate 
then extends to M in the particular case where u° ^ 0. A refined estimate on the 
Ua s is given below in Section 5. Another complement to Theorem 4.1 is that the 
bubbles in this theorem satisfy the ^-decomposition [32]. More precisely, we also 
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have that for any x G M and any ce, 



y(x) = U°(x) + 



Ah,c 



- — ' t ..2 , d g (xj, a ,x) 2 

i= 1 V r'i'OL ' n(n— 2) 



+ Ra(%) 



where the R a ’s, g Hf(M) for all a, are such that R a — > 0 strongly in il 2 (M) 
as a — » +oo. Moreover, if we define the energy of a function u by E(u) = ||ix|| 2 *, 
then 

E(u a ) r = E{u°) r + kK~ n + r a 



for all a, where the r a ’s, r a G M, are such that r a — > 0 as ce — > +oo. In other 
words, the energies split also. 



Atkinson-Peletier [1] and Brezis-Peletier [5] have been concerned with the 
description of the pointwise behaviour of sequences of solutions of equations like 
(1), dealing with radially symmetrical solutions u £ of the semi-critical equations 
A u = on the unit ball of the Euclidean space. More recent developments 

in this specific direction are in Robert [25, 26]. An estimate like in Theorem 4.1 
in the case k = 1, stating that solutions of minimal energy of equations like (1) 
are controlled from above by a standard bubble, appeared then in Han [16] when 
dealing with solutions u £ of the equations A u = u 2 “ 1_£ on bounded open subsets 
of the Euclidean space, in Hebey- Vaugon [19] when dealing with (1) and arbitrary 
Riemannian manifolds, and in Li [20, 21] when dealing with equations like (1) on 
the unit sphere. Improvements, still in the case k = 1, are in Druet [7, 8] and 
Druet-Robert [14]. We refer also to Schoen-Zhang [30]. 



5. Proof of the theorem in the general case 

We only present a very brief sketch of the proof of Theorem 4.1 when we regard 
the theorem under its entire generality. We assume in what follows that 

max u a — ► +oo 

as a — » +oo. If not the case, up to a subsequence, the u Q ’s converge C 2,0 to w°, 
and we do have our asymptotic description. As a remark, the proof of the theorem 
goes through the proof of a slightly stronger result. All what follows is up to a 
subsequence. We let G be the Green function of the operator A g -f /^oo, and let 4> 
be the continuous function on M x M given by 

${x,y) = (n- 2)u> n -id g (x, y) n ~ 2 G(x,y) 



if x t - y, and 



${x,y) = i 
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if x — y, where u; n -i is the volume of the unit ( n — l)-sphere. Then we claim that 
for any converging sequence (x a ) of points in M, and any a, 

k 

u a (x a ) = (l+o(l)^M 0 (x Q ) + y^^(xi,x)+o(l)^S«(a: Q ) 

2=1 

where x is the limit of the x a ’s, X{ is the limits of the a^ a ’s, and B l a is the 
standard bubble with respect to the a^ a ’s and ^ Q ’s given by Theorem 4.1. It is 
easily seen that the theorem and the remarks after this theorem follow from such 
an asymptotic description. The proof of these asymptotics splits into four main 
steps. The first preliminary step consists in getting rescaling invariant estimates. 
Such estimates basically state that there exist k G N*, converging sequences (x^ a ) 
in M, and sequences (//;,«) of positive real numbers converging to 0, i = 1, . . . , fc, 
such that 

R k a (x)^~ l u a (x) < C 

for all x and all a, wheveR^(x) is the minimum over i of the distances from the 
Xi, a ’ s to x. The key idea here is that if such an estimate is false, then we can 
construct another blow-up point. This weak estimate comes with an important 
refinement and complementary informations on the limit of the u a ’s outside blow- 
up points. Then we need to prove that the upper estimate in Theorem 4.1 holds. 
For that purpose, we rearrange the x^ a ’s in families. Inside a family, blow-up 
points are close one to the other. Two families are far one from the other. If the 
y i a ' s are the representative of such families having the largest i = 1, . . . ,p, 
the second step in the proof consists in proving that an upper estimate like in 
Theorem 4.1 holds with respect to the s. We prove that there exist C > 0 and 
R > 0 such that 

u a {x) < (1 + o(l))u°(x) + CJ2 

2=1 

for all a and all x G M\ U p i=1 B Vi where B l a is the standard bubble 
with respect to y^ a and its corresponding The proof of such an estimate 

goes through the establishment of a scale of intermediate estimates, referred to 
as ^-sharp estimates with 0 < 6 < ZL = 1 ^, the weakest of these when 6 = n ^ L 
being like the weak estimate we discussed above. The proof of the upper estimate 
as in the theorem then reduces to the proof that this estimate holds inside the 
B y% a (i?/Xi >Q )’s. This is the third step in the proof. We proceed here by induction. 
We consider subfamilies of blow-up points, and prove that the estimate holds in 
B Vi a (Rni, a )i outside smaller balls of sub-representatives, and so on up to the point 
where we have exhausted all the blow-up points. In each step of this induction 
process, we pass through 6-sharp estimates, and, in some sense, let then 6 — > 0 to 
get the sharp estimate. Once we have proved that the upper estimate of Theorem 
4.1 holds in M, the argument becomes simplier and more conventional. The fourth 
step in the proof consists in proving that the above asymptotics follow from the 
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Green representation formula 

u a (x a ) - u°(x a ) = / G a (x a ,x)(u a (x) 2 -1 - u°(x) r ~ l )dv g 

Jm 

+ G a (x a ,x)(h 00 (x) - h a (x))u°(x)dvg 

JM 

where G a is the Green function for the operator + h a . The different terms 
that are involved in this formula are controlled thanks to the upper estimate we 
have just discussed. The asymptotics follow from rather standard developments. 
As already mentioned, since $ is continuous, the theorem and the remarks after 
the theorem are then easy consequences of the asymptotics. We refer to Druet, 
Hebey and Robert [12] for a complete proof, and also to the next section for the 
special case where k = 1. 



6. Proof of the theorem when the energy is of minimal type 

As a preliminary remark, it is easily checked that the energy of a bubble is a 
constant which is independent of the bubble. More precisely, if B a is a bubble, 
then 

E{B a ) +o(l) 

where o(l) — > 0 as a — > Too, and K n is the sharp constant for the Sobolev 
inequality on the Euclidean space. Returning to equation (E a ), we let (u a ) be a 
sequence of solutions of (E a ). We assume that for any a, 

E(u a ) < tf- ( ”~ 2)/2 • (10) 

The theory in this particular situation is due to Druet and Robert. A com- 
plete proof can be found in Druet and Hebey [10]. Thanks to (10), the u^s are 
bounded in H 2 (M). It follows that there exists u° G Hf(M) such that, up to a 
subsequence, u a — ^ u° in H 2 (M) as a — » Too. Since the energy of the u a 's is less 
than or equal to the minimum energy, see (10), the following holds: 

(PI) either u° ^ 0, and u a — » u° in C 2,e (M) as a — » +oo, 

(P2) or u° = 0, and the u a ’s blow-up with one bubble. 

The general idea here is that we do not have enough energy so that we can face a 
multiple bubbles situation, neither the situation where we would have a nonzero 
limit u° and a bubble. Now the statement we want to prove is the following. We 
assume that u° = 0, if not we do have our asymptotic description, and we assume 
that (10) holds. Then we want to prove that there exists some constant C > 1, 
that there exists a converging sequence (x a ) of points in M, and that there exists 
a sequence (/i a ) of real numbers, fi a > 0 and jj, a — » 0 as a — > +oo, such that, up 
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to a subsequence, 



1 

C 



l~^a 






dg ( I Q ,x ) 2 

n(n — 2) 



n — 2 
2 



< u Q (x) < C 



f^a 



X a + 



d g (x a ,x) 2 
n(n — 2) 



n-2 

2 



( 11 ) 



for all x, and all a. This is exactly Theorem 4.1 in this particular situation where 
the energy is assumed to be of minimal type. In such a particular situation, we 
can define explicitely the x Q ’s and the /i a ’s. More precisely, the x a ’s and /i^s are 
here given by 



l — — 

u a (x a ) = maxii Q = fia. 2 
M 



( 12 ) 



In other words, x a is a point where u a is maximum, and fi a is, up to some power, 
the inverse of the maximum of u a . As when discussing the general case, we need 
first to prove the upper estimate in our statement, namely in (11). A simple claim 
then, thanks to (12), is that this upper estimate is equivalent to the following 
estimate: there exists C > 1 such that, up to a subsequence, 



_ n — 2 

d g (x a ,x) n ~ 2 ii a 2 U a (x) < C 



(13) 



for all x, and all a. In order to prove (13) we proceed in three steps. All what 
follows is up to a subsequence. As above, the first step consists in proving rescal- 
ing invariant estimates. To some extent, these estimates state that the are 
controlled by the square root of a bubble. They can be stated as follows: 



Step 1 (Rescaling invariant estimates): There exists C > 0 such that for any a 
and any x, 

c( g (x Q ,x)“ u a (x) < C . (14) 

Moreover, for any R > 0, for any a , and for any x E M\B Xa (R/i a ), 

d g (x a ,x) n ^ 1 u a (x) < e R , a (15) 

where £r, q is such that lim lim = 0. 

/?— > + oo a— > + oo 

We have here two estimates. The second estimate is a refinement of the first one. 
The general idea of the proof of (14) and (15) is that if one of these two estimates 
is false, then we must have another bubble in the decomposition of the u a ’s. This 
is not possible since we assumed that the energy of the u a ’s is of minimal type, 
so that we just have one single bubble in the decomposition. In Step 2, we prove 
that the upper estimate in our statement (11) holds up to some £, 6 > 0. In other 
words, we prove something like an e-sharp upper estimate in (11). 

Step 2 (e-sharp upper estimate): For any e > 0, there exists C £ > 0,such that 

dg(x a ,x) n ~ 2 ~ E Ha~ +£ U a (x) < C £ (16) 

for all x, and all a. 



As a remark, (14) gives that (16) holds with e = (n — 2)/2. Another remark 
is that (16) becomes stronger and stronger as e decreases, e — > 0. The proof of 
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this estimate makes use of the existence of a Green function for the second order 
operator in the left-hand side of equation (1). What we actually prove, thanks to 
maximum principles, is that the s are bounded from above by a (1 — 6)-power of 
this Green function. This implies (16). In Step 3, we prove that the upper estimate 
in our statement is also true. In other words, we prove that we can take e = 0 in 
the e-sharp upper estimate (16). 

Step 3 (Upper estimate): There exists C > 0 such that 

dg(x a , x) n ~ 2 fia~ U a (x) < C (17) 



for all x, and all a. 

As already mentioned, (IT) is exactly the upper estimate of our statement (11), 
namely that the s are bounded from above by some constant C times the 
standard bubble defined by the x a ’s and /i a ’s. The proof of this estimate goes 
through the Green representation formula, the s being controlled thanks to the 
6-sharp upper estimate (16), 6 > 0 small. With such an estimate we have half of 
our statement. Summarizing Steps 1-3, we prove first that we do have rescaling 
invariant estimates. This corresponds to the case where e — (n — 2)/2 in (16). 
Then we prove that the 6-sharp upper estimates (16) hold when 0 < 6 < (n — 2)/2. 
Then, at last, we prove that our upper estimate (17) holds. This corresponds to 
the case where 6 = 0 in (16). Step 2 uses Step 1. Step 3 uses Step 2. As a remark, 
we assumed here that the h a ’s converge in C 0,e to some limit function hoo. The 
rescaling invariant estimates actually hold if we only assume that the negative 
part of h a is small in some L p -norm. The 6-sharp upper estimate and the upper 
estimate hold if we only assume that the h a ’ s are greater than or equal to some 
function V for which the operator + V is coercive. 

Now that we have the upper estimate in (11), we can discuss the last step in 
the proof. We let B a be the standard bubble in (11), and G be the Green function 
of the operator A g + h^. We let also xo be the geometrical blow-up point, defined 
as the limit of the x a ’s as a — > +oo. At last, as in the preceding section, we let 
: M x M -> M be given by $(x,y) = (n - 2)o; n _id^(x, y) n_2 G(x, y) if x ^ y, 
and <h(x,y) = 1 if x = y, where cj n _i is the volume of the unit (n — l)-sphere. A 
simple claim, which follows from the construction of the Green function, is that $ 
is continuous on M x M. Step 4 in the proof can then be stated as follows: 



Step 4 (Asymptotic behaviour): for any converging sequence (x a ) of points in M, 



a 



lim 



Ug{X a ) 

+oo B a (x a ) 



= $(x 0 ,xo) 



where xq is the limit of the x a ’s as a — > -foo. 



We prove Step 4 by using once again the Green representation formula, the u Q ’ s 
being controlled thanks to the upper estimate of Step 3. As already mentioned, 
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our conclusion easily follows from Step 4. Noting that <f> is positive, and since 
is continuous, we get with Step 4 that there exists C > 1 such that 

-jjBa(x) < U a (x) < CB a (x ) 

for all x, and all a. This proves (11), and thus Theorem 4.1 in this particular 
situation where the energy is assumed to be of minimal type. Moreover, since 
<h = 1 on the diagonal, and since <f> is continuous, we also get with Step 4 that for 
any e > 0, there exists 5 e > 0 such that 

j^B a (x) < u a (x) < (1 + e)B a (x) 

for all x E B Xo (5 £ ), and for all a. This proves, in this particular situation where the 
energy is assumed to be of minimal type, the remark after Theorem 4.1 concerning 
the sharpness of the constant. 
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Application of Scans and 
Fractional Power Integrands 

Thierry De Pauw and Robert Hardt 



Abstract. In this note we describe the notion of a rectifiable scan and con- 
sider some applications [DH1], [DH2] to Plateau-type minimization problems. 
“Scans” were first introduced in the work [HR1] of Tristan Riviere and the 
second author to adequately describe certain bubbling phenomena. There, the 
behaviour of certain W 1,3 weakly convergent sequences of smooth maps from 
four-dimensional domains into S 2 led to the consideration of a necessarily in- 
finite mass generalization of a rectifiable current. The definition of a scan is 
motivated by the fact that a rectifiable current can be expressed in terms of 
its lower-dimensional slices by oriented affine subspaces. By integral geometry, 
the slicing function for the rectifiable current is a mass integrable function of 
the subspaces. With a scan one considers more general such functions that 
are not necessarily mass integrable. 



1. Rectifiable currents and the plateau problem 

An m-dimensional rectifiable set R in R n is a subset of some countable union 
U iZoMi where Mi, M 2 ,... are m-dimensional C 1 submanifolds and Mo has m- 
dimensional HausdorfF measure W m (Mo) = 0. At TV 71 almost every point x E i?, 
R has an approximate tangent space. [S], Section 3. An m-dimensional (integer- 
multiplicity) rectifiable current T in R n is given by a bounded m-dimensional 
Borel measurable rectifiable concentration set Rt together with an 7^ m -integrable 
density function 0t : Rt — * {1,2,...} and an 7-f m -measurable orientation T : 
R T ->• Am R " so that at H m almost every x E Rt , T(x) is the wedge product of 
vectors from an orthonormal basis of the approximate tangent space of Rt at x. 
See [FI] ,4. 1.24 or [S], Section 27. Thus the action of the current T on a differential 
m-form <j> E © m (R n ) is given by the integration 

T(4>) = [ (f(x),<j>(x))e T (x)dH m x . 

J Rt 

The mass of T is then simply M(T) = J R 6t{x) dhC^x. For m > 1, boundary of 
T is the m — 1-dimensional current defined by the formula dT(ip) = T(dip) for 
ijj E V ni ~ -HR"). A rectifiable current generalizes an oriented submanifold M. We 
sometimes use the abbreviated notation [M], in case the orientation is known, 
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for the corresponding multiplicity one rectifiable current. So if M is an oriented 
manifold with boundary, Stokes Theorem becomes d\M\ = [<9MJ. Even though 
H m (RT ) < oo, one should be aware that the support of the current, sptT, which 
may be much larger than Rt , might even be n-dimensional for m > 1 . On the other 
hand a zero-dimensional rectifiable current is simply a finite integral combination 
of point masses. Let TZ m denote the group of ra-dimensional rectifiable currents 
in R n . In 1960, H. Federer and W. Fleming obtained the following fundamental 
existence theorem: 

1.1 Theorem. [FF] Given any To G 7 Zm with dTo G the family of currents 

{T G IZm • dT = dTo} contains a rectifiable current of least mass. 

This theorem is valid for all m > 1 in any R n as well as in any compact Riemannian 
manifold (provided the admissible family is nonempty). There one also has, in 
any homology class, a rectifiable current that minimizes mass. Among general 
currents the existence of a mass minimizers is an easy consequence of the Banach- 
Alaoglu theorem, but what is important in [FF] is the rectifiability, which should be 
understood as an initial regularity for minimizers. The complete interior regularity, 
of such rectifiable mass-minimizers, i.e., that sptT \ spt<9T is an embedded real 
analytic submanifold, was established in the sixties for 1 < m = n— 1 < 6 by works 
of Fleming [FI], De Giorgi [D], Almgren [Al], Triscari [Tr], and Simons [Ss] . De 
Giorgi’s work showed that a seven-dimensional mass-minimizer in R 8 would have 
at most isolated interior singularities, and, in fact in 1970, Bombieri, De Giorgi 
[D] , and Giusti [BDG] established the mass-minimality of the specific example 

Q = d[{{x,y) e R 4 X R 4 : \x\ < |y|}J L_ Bf , 

which has an isolated singularity at (0, 0). Then the cartesian product with a cube 
Q x [—1, 1] J is mass minimizing in R 8+J \ and so the following result of H. Federer 
gives the optimal estimate of the Hausdorff dimension of the interior singular set 
of a codimension one minimizer: 

1.2 Theorem. [F2] For any m-dimensional mass-minimizing rectifiable current T 
in R m+1 and e > 0 



^m- 7 +e ( s^g ( S pt T \ spt dT)) = 0. 

The complete boundary regularity of T was established by Hardt and Simon [HS] 
near any point where the given dTo is a smooth m — 1-dimensional oriented em- 
bedded submanifold of R m+1 . In 1984, F.J. Almgren completed a massive work 
treating the higher codimension interior partial regularity: 

1.3 Theorem. [A4] For any m-dimensional mass-minimizing rectifiable current T 
in R n and e > 0 

W m — 2 +e ( Sing ( sp tT\ sptST)) = 0, 

and S. Chang [C] showed that interior singularities are at most isolated points for 
m — 2. Note that [FF] already contained the singular minimizing example of the 
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sum of two oriented totally orthogonal disks in R 2 x R 2 . 



[B? x {0}] + [{0} x B 2 ] . 



2. Size and fractional powers of the density 

In R 3 , two-dimensional mass- minimizing rectifiable currents have no interior sin- 
gularities and provide a nice model for some but not all “soap films” . General soap 
films may have interior singular curves which simultaneously border three surfaces 
meeting at equal angles. To use currents in a better model for soap films, Almgren 
[A3] introduced the notion of size for a rectifiable current: 

Size (T) = n m (R T ) • 



Thus one ignores the density function in computing size. To understand size versus 
mass minimization, consider the one-dimensional rectifiable current in the plane 
consisting of 2 parallel similarly oriented intervals 



T n = K-- 



2 ’ 



A (I Al + „_i _ A ,1 

2 M 2’ 2 ^ ^ ^ 2’ 2 M 2’ 




Then To is mass-minimizing among all rectifiable currents having boundary equal 
to <9 Tq. However, the size minimizer is sum of five oriented intervals, 



Ti = I(-|,^),(0,0)] + (0,0)] 

+ 2[(0,0), (3,0)] + 1(3,0), (-, T)J -|- [(3,0), (-, T)] , 
one of which has multiplicity 2. Note that 

M(T 0 ) = 2-4 < 4-1+2-3 = M(Ti) and Size(T 0 ) = 8 > 4-1+3 = Size(Ti). 



Similarly in dimension 2 one may consider the sum To of two close coaxial, parallel 
and similarly oriented, disks. Then To is mass-minimizing, but the size-minimizer 
T\ with boundary dTo contains a single multiplicity 2 disk in the middle and 
the set spt T\ models a soap film with an interior singular curve. A fundamental 
problem with size-minimization is the lack of a general existence theorem. One 
is faced with the possibility of a size-minimizing sequence of rectifiable currents 
having unbounded masses and failing to have subsequences convergent as currents. 
This is what happens in an example of F. Morgan. 



2.1 Example [M]. For a fixed 1 < (3 < 2 consider the following countable sum of 
vertical oriented intervals in the plane: 

00 1 1 11 

Then T 0 G because M(T 0 ) = Yl'jLi fv < 00 • One may check that To is mass- 
minimizing by an easy callibration argument. However, one sees, as in the previous 
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section, that one may decrease the size, at the expense of increasing the mass, by 
replacing the oriented interval [(1,-1), (1,1)] by the sum of three intervals 

[(i,-i), (l,i)] + i(l,-^),(^)] + l(^), (i,i)]. 

The new current has the multiplicity 2 interval 2|(|, — ^), (|, ~p)j which we may 
then replace by the sum 

2[(l,-^),(l,-p)l + 2[(l,-p),(l,p)J + 2[(1,T),(I,-^)]. 

Continuing we obtain a size-minimizing sequence whose mass approaches infinity. 
These do not converge as currents and the resulting formal countable sum of 
oriented intervals with integer multiplicities is not a current. We will see that 
this can be understood as a “scan”. There are some positive results concerning 
size-minimization. 

2.2 Theorem. [M] If spt (dTo) is a smooth m — 1- dimensional submanifold of 
R m+1 that lies on the boundary of a smooth compact convex body, then there 
exists a size- minimizing rectifiable current T with dT = dTo- 

The idea here is that one can modify a size-minimizing sequence using decomposi- 
tions into oriented boundaries of sets which extend all the way to the boundary T 
of the convex body. Since T \ spt dTo contains only finitely many components, one 
thus obtains a bound Ti™ almost everywhere on the densities in the sequence. So 
the masses are bounded, and one has convergence as currents. The lower semicon- 
tinuity of size under this condition was established in [A3]. In [DH1] we obtain, for 
a general codimension one rectifiable current Tq , (without the convex hull prop- 
erty of spt (dTo)) a weaker result concerning the existence of a minimizing set. An 
ra-dimensional set S C R n is minimizing [A2] relative to a compact set K if 

n m [f(s)\ > n^is) 

for any Lipschitz map / : R n — ► R n with {x : f(x) / x} C R n \ K. 

2.3 Theorem. [DH1] If, in R m+1 , To G Hm and 7Y m ( spt <9T 0 )) — 0, then there 
exists a minimizing set S with respect to spt dTo having spt dTo C S. 

Here is a brief outline of our construction. First we penalize the lack of compactness 
by choosing, for any 0 < e < 1, a rectifiable current T e minimizing Size(T) + 
eM(T) among rectifiable currents T with dT e = dTo. Second we observe that the 
(renormalized) measure 

He = ( H m L_ i?T e )(l + 

defines a (real) rectifiable varifold, stationary [Al] in R m+1 \ spt dTo ■ That is, the 
first variation of the total mass of the varifold fi e vanishes for any deformation 
which fixes spt<9To- Since the total mass of each \i e is bounded by Size(T 0 ) + 
M(To), we may by [Al], choose a sequence pi ej weakly convergent to a rectifiable 
varifold fi. We show that the m density of this varifold is one at H m almost 
every point of spt /i \ sptSTo- Thus /a = H 171 l_ S for some rectifiable set S. 
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Third, by modifying some arguments of Ambrosio, Fusco and Hutchinson [AKH], 
which proved the minimality of a limit of codimension one minimizing sets, we 
verify that S is the desired minimizing set. It is an interesting question to find 
conditions that will guarantee that the currents T ej converge. While Morgan’s 
example above shows that this will not happen in general, the condition that 
spt <9To be a smooth submanifold may be sufficient. Even for the “soap-film case” of 
two-dimensional size-minimizers in R 3 , lack of a priori knowledge of the boundary 
behaviour prevents progress. In the next result, one avoids the boundary behaviour 
problem. 

2.4 Theorem. [M], [DH1] In a compact Riemannian 3-manifold, any two-dimen- 
sional homology class contains a size-minimizing rectifiable current. 

J. Taylor [T] classified the local interior structure of the minimizing set S. Up 
to a C 1,a diffeomorphism of space, the neighborhood of an interior point is ei- 
ther a plane, three half-planes meeting at equal angles along a line, or 6 planar 
sectors meeting at equal dihedral angles at a point. This local structure and the 
compactness of S imply that there exists a global Lipschitz retraction p of some 
open neighborhood of S onto S. Since sptT ej = spt p ej converge in Hausdorff 
metric to S, we find that, for some fixed, sufficiently large j, the retracted current 
p#T ej is the desired size-minimizing homology representative. To use this retrac- 
tion argument for the original Plateau boundary problem, one would need more 
information about the boundary behaviour of minimizing sets or (approximately) 
size-minimizing currents. One can also consider a free boundary or obstacle prob- 
lem where the minimization occurs among rectifiable currents constrained to have 
their boundaries lie on a given hypersurface. Here, for two dimensions in three 
space, the local structure is easy to classify. Again up to aC 1,Q diffeomorphism of 
space, the neighborhood of a free boundary point is either a half-space perpendic- 
ular to the given hypersurface or else 3 quadrants at equal angles to each other, 
meeting the hypersurface orthogonally. Thus one obtains 

2.5 Theorem. [DH2] For any smoothly bounded region ft C R 3 and any To G 1Z 2 
with spt To C R 3 \ ft and spt d To C dft, the relative homology class 

{T G IZ 2 : sptT C R 3 \ ft, spt3T 0 C dft , and 

spt (T — T 0 — dS ) C dft for some S G 7^3 } 
contains a rectifiable current of least size. 

Note that one can imagine ft as the interior of a thick three-dimensional wire 
whose boundary surface is supporting the boundary of a soap film. Unfortunately 
this last existence theorem provides no obvious uniform mass bounds on the size- 
minimizing currents for a sequence of such thick wires fti shrinking to a one- 
dimensional smooth curve. There are many other functionals that may also give 
rise to some minimizing sequences which have unbounded mass and which have 
no subsequences convergent as currents. Perhaps the simplest are given by “frac- 
tional” powers of the density. 
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2.6 Definition. For any a G [0, 1] and T G 7 £ m , the a-mass of T equals 
M a (T) = [ [0 T (x)] a dH m x . 

J Rt 

The case a = 0 corresponds to the Size(T) = H m (RT) and the case a = 1 
corresponds to the ordinary mass M(T). The special property of the range 0 < 
a < 1 is that the functional M a is weakly lower semicontinuous on mass bounded 
sequences. In case m = 0 one is only dealing with “collisions of atoms” and the 
result is elementary, based on the inequality 

\i+j\ a < |i| a + \j\ a for any integers i and j , 

which is only valid for a G [0, 1]. For larger m one may reduce to the case m — 0 
by slicing which we now review. 



3. Slicing 

For any bounded Borel m form Q, on R n , let T I Q denote the zero-dimensional 

current given by (T l_ £i)(ip) = Ti^VL) for ^ G Co°(R n ). From [FI], 4.3 we recall 
that if / : R n — > R m is Lipschitz and T G 7 Z m , then, for a.e. y G R m the slice 
(T, /, y) e where 

(' T,f,y ) = limT l_ with Cl m — dx 1 A • ■ • A dx m . 

N ' rjO v c^ m r m 7 

In fact, for a.e. y G R m , Rt n/ _1 {z/}isa finite set of points x where the approx- 
imate tangent space of Rt exists with f\Rr being approximately differentiable of 
rank matx, and the slice is given by the formula 

( T,f,y ) = a(x)0 T (x)S x 

xeRr^f~ 1 {y} 

where a(x) = sgn ( /\ m Df(x)T(x ), fi m ). One also has the integral formulas 

[ (T,f,y)dy = T L f*Q m , 

J R m 

[ M a (TJ,y)dy = M a [T l_ f*n m ] • 

J n™ 

Recall also the compact space V of orthogonal projections of R n onto R m , which 
has an invariant probability measure induced by the transitive action of O(n). In 
particular, for each increasing function A : {1, . . . , m} — > {1, . . . , n} we have the 
coordinate projection 

PX : R n -* R m , p\{xi,...,x n ) = Oz A(1 ),...,:r A ( m) ) . 

By writing an m form (j) G V (R n ) in coordinates 

<f> = ^ 4>\dx A = Y^xpftlm 

A A 
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with each G Co°(R n ), we see how the current T is expressed in terms of its 
zero-dimensional coordinate slices 

T{4>) = J2 T = Z( T| -pfn m )(4>x) 

A A 

= V] [ {T, p\, y){4>\) dy ■ 

Composing with rotations and integrating over O(n), we find a formula 

T(<j>) = (volO(n)) -1 [ [ { T - l P\° 9 ,y){ < t>-{P\° 9 )*^m)dydg 

JO(n) A J R m 

(volO(n)) -1 f [ (T,p,y)(<t> -p*£l m )dydp . 

Jv J R™ 

One also has the integral geometric relation 

M a (T) = (3(m,n) [ [ M a (T,p,y) dpdy . 

JV j R m 

This motivates us to define various classes of scans as measurable functions 
T : V x R m -f { zero — dimensional currents } . 

In particular, a rectifiable scan is a measurable function 

T : Px R m — » IZo 

corresponding to some rectifiable set Rt , some 7f m -measurable 6t • Rt — >• 
{1,2,...} and some W m -measurable orientation T of Rt so that, for almost all 
(M)^x R™, 

T(P,y) = <t(x)0 t (x)6 x 

x£R T C\p- x {y} 

where a {pc) = sgn ( f\ rn pT{x), fi m ). One defines 

M a (T) = P(m, n) f f M a (T(p,y)) dpdy , 

JV J R m 

and sees that the scan T corresponds to a rectifiable current if and only if Mi (T) < 
oo. For a rectifiable current T G 7£ m , an elementary Fourier transform argument 
shows that 

dT = 0 if and only if (T,p, p)( 1) = 0 for almost all (p,y) G?x R m . 

Thus, for 2 rectifiable scans 5, T we say 

dS = dT if and only if (S — T)(p,y)( 1) = 0 for almost all (p, y) . 

3.1 Theorem. [DH1] Suppose 0 < a < 1, To G 7 Z m , & n d W m (spt<9To) = 0. Then 
there exists a rectifiable scan T with dT — dTo and 

M a (T) = inf{M a (T) : T G 7£ m , dT = dT 0 } . 
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The proof involves getting the convergence as scans of a subsequence of an M a 
minimizing sequence. We must work with a convergence that is weaker than the 
weak convergence of currents. To see this, consider the following 

3.2 Example. Suppose Sj — j 2 <9[Bij] in R 2 . Thus Sj is an oriented circle of 
radius i with multiplicity j 2 . Then, in the weak topology of currents, 

Tj ->■ <9[<5( 0) o)ei A e 2 ] ^ 0 , 

but, for 0 < a < 

0 HT 

Ma(T,-) = j 2a -(j) - 0. 

An appropriate topology is the a flat distance [F12], [W] on R m 

X a {T u T 2 ) = inf{M a (R) + M a (S) : T 1 -T 2 = R + dS, R E Km, S E ft m+1 } • 

Recall now from [AK] that a measurable function / from R m into a general metric 
space X has finite total variation if 0 o / is BV for all Lipschitz <f : X — > R and 

||£>/||(R m ) = sup{/ TO o/)| : (ft : X —> R, Lip</»<1} < oo . 

R. Jerrard [JS] observed that the slice of a normal current was of metric bounded 
variation with respect to the flat norm. The next result is the analogue for the T a 
distance. 

3.3 Theorem. [DH1] IfT E Rm, dT E M a (T) + M a (9T) < oo, andp E V , 

then 

(K 0 ,X a ), f(y) = ( T,p,y ) 

has finite total variation 

||D/||(R m ) < m [M a (T) + M a (<9T) ] 

Proof. In case m — 1 the equation, for almost all 5 < t, 

(T,p,s) - ( T,p,t ) = (dT) L_ p^[s,t] - d(T l_ p _1 [s, t}) 
implies that 

XF a ((T,p,s),(T,p,t » < M a (T l_ p -1 [s,£]) + M a ((dT) Lp-'lM]) • 
Summing over almost all partitions of R then gives that the essential variation 
essvar (/) < M a (T) + M a (<9T) , 

which implies the case m = 1. For m = 2, 3, . . . , we use the formula [FI], 4.5.9.27 

[ \m 

J R™ 

m r. 

< T / essva,Tip(y 1 ,...,y i - 1 ,-,yi +1 ,...,y m )dy 1 ...dyi- 1 dy i+ i...dy m 

i=i jRm ~ 1 
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with = (j)o f to get the desired bound. This estimate gives the desired variation 
in y. To get a similar variation bound in (p, y) one may use, for almost all (p,y), 
the formula 

(T,p,y) = U#(Tx[V],P,(y,p)) 

where U(x,q) = x and P(x,q) = {q{x),q) for (x, q) G R n x V. From the usual 
BV compactness for real- valued functions it is not difficult to derive the following 
metric-space-valued version: 

3.4 Theorem. [DH1] Suppose N is a Riemannian manifold, Y is a separable 
metric space, M : Y — > R + is lower semicontinuous, and M~ 1 [0,R] is sequentially 
compact for all R > 0. If fj : N — > Y is measurable with 

A = sup||£>/j||(iV) + [ M(fj(x )) < oo, 

3 JN 

then some subsequence fj> converges pointwise a.e. to a measurable f : N — > Y 
with 

\\Df\m + [ M{f(x)) < A . 

JN 

To prove the existence of a suitable scan T for Theorem 3.1, we can now apply 
Theorem 3.3, the remark after, and Theorem 3.4 with 

N = Pxn m ,Y = (K 0 ,P a ), M = M«, fj(p,y) = (T^y). 

The lower semicontinuity follows from Fatou’s lemma and the lower semicontinuity 
of M a on (T^o j^a)- To verify the rectifiability of the limiting scan, we obtain a 
rectifiable varifold as in the proof of Theorem 2.3 and then show that our scan 
is necessarily concentrated on an 7f m -finite rectifiable concentration set of this 
varifold. Almost every projection is transverse a.e. to this rectifiable set and one 
works with the convergence of the slices in these directions to eventually get the 
desired consistently defined multiplicity and orientation T. In [DH2] we obtain 
a general compactness theorem in the class of rectifiable scans having rectifiable 
boundaries and having M a + M a <9 uniformly bounded. This gives an existence 
theory for various Plateau problems. We also obtain the optimal interior partial 
regularity estimate for M a minimizers: 

3.5 Theorem. [DH2] For any m- dimensional M a minimizing rectifiable scan T in 
R n and e > 0 

T^m-i+e ( s^g ( sp t T \ spt dT)) = 0. 

Moreover, T \_ K has finite mass for every compact K C R n \ spt dT. 

The proof uses the following (roughly stated) 

3.6 Lemma. IfO < a < 1, then, near any point a of spt T\ spt dT where T has a 
multiplicity v tangent plane, the scan T is mostly a multiplicity v graph of a single 
function. 
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In shrinking cylinders about such a point, the v excess [FI], 5.3, as well as the 
ordinary excess, approaches zero. This is in sharp contrast to the case a = 1 of 
mass-minimizing currents. For example, the complex cusp 

{(z, w) E C x C : w 2 = z 3 } 

supports a multiplicity one mass-minimizing rectifiable current with a multiplicity 
two tangent plane at (0,0). The proof of Lemma 3.6 is based on a multi-valued 
graphical approximation followed by a squashed comparison current. 



4. Another fractional integrand 

In the works [HR1], [HR2] treating various energy-bounded sequences of Sobolev 
mappings, one encounters rectifiable currents with bounds on the integral of frac- 
tional powers of the (ordinary) mass of slice. In our notations, the analogous 
situation is to consider for T G 7 £ m and 0 < a < 1, the integral 

M a (T) = (3(m,n ) f f (M {T,p,y)) a dpdy . 

JV J R m 

Note the inequalities 

M a (T) < M a (T), M(T)<M!(T)<M(T) . 

4.1 Example. Consider in R 2 the concentric, multiplicity one circles 

j 

T i = J2 d l B i/il efti- 

Then, as j — » oo, 

M q (Tj ) = M(Tj) — > oo , but supM a (Tj) < oo . 

j 

Here the limit is a rectifiable scan with M a finite, but with a concentration set of 
infinite Hausdorff measure. As in [HR1], [HR2], one may again work with (7£o> •?■<*) 
to obtain the lower semicontinuity of M a . One can obtain the existence of M a min- 
imizing scans by using a weak- type bound in the Lorentz space L« ,oc . See [HR1], 
Section 9. However the partial regularity or even rectifiability of M a minimizers 
is unknown. See [HR1], Section 8.2. 

5. Related problems 

In the study of M a minimizers, one may replace the power function 6 a by a smooth 
concave unbounded increasing function H(6) with H( 0) = 0 and 1) = 1. Use 
of scans accommodates as well treatment of the case when H = H(Q,x) is also 
allowed to depend smoothly on the space variable x. By using H(9) as an alternate 
norm on the group of integers, a rectifiable minimizer may also be found in a 
generalized class of flat chains following the works of Fleming [F2] and White [W] . 
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The close relation between rectifiability and slicing explained in [AK] and [W] 
was an important motivation for the definition of scans in [HR1] and [DH1]. One- 
dimensional flat M a minimizers are also applied to describe transport paths in 
[XI], [X2]. The M a functional reflects the efficiency of combining paths in various 
distribution systems, such as mail delivery, the circulatory system, etc. Here one 
obtains, for any two probability measures /io, /jl\ in R n , a flat 1 chain T which 
minimizes M a subject to the constraint dT = /ii — /io as zero-dimensional currents. 
In [XI], the transport path T has positive real density function Or with values in 
(0, 1] because M(/ii) = 1 = M(/i 2 ). So one here has the inequality 

M a (T) > M(T) , 

which is just the opposite of the inequality that we had in our study [DH1], [DH2] 
of integer- multiplicity rectifiable currents. In contrast to [DH1], [DH2], M a min- 
imization for transport paths always gives finite mass currents. In [X2], Q. Xia 
proves the precise local interior regularity of an M a minimizing path: that, in 
R n \ ( spt fi\ U spt /X 2 ), T is locally a finite collection of oriented intervals with mul- 
tiplicities. Finally we have begun work in [DH3] on carrying over various results 
of geometric measure theory to rectifiable scans. There we introduce the .notion of 
a rectifiable scan in a metric space X, and give some results generalizing the work 
of L.Ambrosio and B. Kirchheim [AK] on currents in metric spaces. The idea is 
that an m-dimensional rectifiable scan in X is a measurable function 

T : Lip (X, R m ) x R m -> U 0 (X) 

which admits a representation in terms of an m-dimensional rectifiable [AK] set 
Rt C X, and integer density function and an orientation T of R^-. (The notion 
of orientation requires some effort to describe.) For a Lipschitz map g : X — » Y of 
metric spaces and m-dimensional rectifiable scan T on X, the push- forward g#T 
defined by 

(a#T)(f,y) = g#(T(f°g,y)) 

for almost all (/, y) £ Lip(X, R m ) x R m , is a rectifiable scan on Y. An m — 1- 
dimensional scan S is the boundary , 9T, of T if, for all points a £ X, 

S{g,z) = limT((fl r , dist(-,a)),(r,z)) 

r jO 

for almost all (g,z) £ Lip(X, R m_1 ) x R m_1 . One may again impose bounds on 
M a and M a <9 along with suitable topological bounds on supports to obtain scan 
compactness theorems. 
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Bubbling of Almost-harmonic Maps between 
2-spheres at Points of Zero Energy Density 

Peter Topping 



Abstract. We show that bubbling of almost-harmonic maps between 2-spheres 
has very different behaviour depending on whether or not bubbles develop at 
points in the domain at which the energy density of the body map is zero. 
We also see that this translates into different behaviour for the harmonic map 
flow. In [11] we obtained results, assuming nonzero bubble point density for 
certain bubbles, forcing the harmonic map flow to converge uniformly and 
exponentially to its limit. This involved proving a type of nondegeneracy for 
the harmonic map energy (a ‘quantization’ estimate’) and an estimate on 
certain bubble scales (a ‘repulsion’ estimate). Here we show that without the 
nonzero bubble point density hypothesis, both the quantization and repulsion 
estimates fail, and we construct a flow in which the convergence is no longer 
exponentially fast. 



1. Introduction 

In this article, we study the harmonic map energy functional 

E ( u ) -= \ ( I Vm I 2 

^ Js 2 

for maps u : S 2 — » S 2 c -> M 3 . We are interested in the behaviour of E near 
its critical points - known as harmonic maps - and near more general objects - 
known as bubble trees - which can be thought of as critical points at infinity. Such 
information is related to the behaviour of the L 2 -gradient flow for E ~ known as 
the harmonic map flow - which is a solution of 

^ = T(u(t)), (1.1) 

where the tension T{u) G T(u*(TS 2 )) is defined to be the negation of the L 2 - 
gradient of E. (See [11] or [4] for more details in this case.) The flow, when smooth, 
dissipates energy according to 



jE(u{t)) = -\\T{u{t))\\l H s^ 
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which integrates to the useful identity 

/ oo 

\\ T ( u (t))\\LHs*) dt = E ( u ( s )) - E i u (t)) < oo. (1.2) 

Harmonic maps between 2-spheres are well known to be precisely the rational 
maps and their complex conjugates (see [4, (11.5)]). In particular, being conformal 
maps from a surface, their energy is precisely the area of their image, and thus 

E(u) = 47r| deg('u)! £ 47rZ, 

for any harmonic u. Since the harmonic maps occur in families, the energy E is 
degenerate near each critical point. However, by combining work of Gulliver- White 
[5] and L. Simon [7], [8, Lecture 3, Section 2] it is possible to prove a substitute 
for nondegeneracy, that 

0 < E(u) - E(h) < C\\T(u)\\l HS2) , (1.3) 

for u sufficiently close 1 to a harmonic map h : S 2 — > S 2 . Estimate (1.3) can be 
used to analyse the asymptotics of the harmonic map heat flow when it converges 
smoothly (without bubbling as described below) at infinite time. In particular, for 
a flow converging smoothly to a limit u as t — > oo, 

§j. iE(u(t)) - E(u x )} = -||T(u(t))||£ 2(52) < -i [E(u(t)) - E ( Uoo )] , 

for sufficiently large t , and we see that E(u(t)) — E(u oo) must decay to zero expo- 
nentially fast. Moreover, (1.3) allows us to calculate 

[E(u(t)) - E( Uoo fi = i [E(u(t)) - £(««,)]-* \\T(u(t))\\l HS2) 

~ 2v^H' ? "( U ( i ))^ L2 ( 52 )’ 

and thus, for s £ [0, oo) (and a new constant C independent of s) 

poo 

/ ||T(«(0)||l=(^)^<C[£( W (s))- j B( Woo )]5. (1.4) 

To begin with, this implies the finiteness of the left-hand side of (1.4) (which is 
simply false for the harmonic map flow in general, despite (1.2)). In addition, 
the exponential energy decay then forces the left-hand side of (1.4) to decay to 
zero exponentially fast. By returning to the equation (1.1) we then must have 
exponential convergence of u(t) to u Q 0 as t — > oo in L 2 (S 2 ) and thus in C k - for 
any k - by interpolation. (Even the convergence of u(t) to u Q 0 in L 2 (S 2 ) uniformly 
as t — > oo is false in general for the harmonic map flow - see [10].) Unfortunately, 
in general the harmonic map flow does not converge smoothly at infinite time (see 
[1] or the examples below). The best we can hope for is ‘bubbling’ convergence as 
the flow approaches a ‘bubble tree’ which can be thought of as a critical point at 



1 for example, close in C 3 
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infinity for E. We first describe bubbling in the context of sequences of almost- 
harmonic maps. The following theorem is a combination of the seminal work of 
Struwe [9] and the improvements of Qing [6], Ding-Tian [3], and Wang [12]. See 
[11, Theorem 1.1] for a more sophisticated version. 



Theorem 1.1. Suppose that u n : S 2 -> S 2 ^ M 3 is a sequence of smooth maps 
which satisfy E(u n ) < M for some constant M, and all n G N, and T(u n ) — > 0 in 
L 2 (S 2 ) as n — > oo. Then we may pass to a subsequence in n, and find a harmonic 
map Uoo : S 2 — > S 2 , and a set {x 1 , . . . ,x m } C S 2 such that 

(a) u n — ^ Uoo weakly in W 1,2 (S 2 ), 

(b) u n -► Uoo strongly in Wf^{S 2 \{x l , . . • ,x m }). 

Moreover, for each x j , if we precompose each u n and Uoo with an inverse stere- 
ographic projection sending 0 G M 2 to x j G S 2 (and continue to denote these 
compositions by u n and Uoo respectively) then for i G {1, . . . , k} (for some k G N 
dependent on j ) there exist sequences a l n — > 0 G M 2 and X l n { 0 as n — » oo, and 
nonconstant harmonic maps to 1 : S 2 -> S 2 (which we precompose with the same 
inverse stereographic projection to view them also as maps M 2 U {oo} — ► S 2 ) such 



that: 

(i) 



x i K x^xi 



as n — > oo, for each unequal i,j G {1, . . . , k}. 



k 

(ii) lim lim E(u n , D^) = E(w l ). 

Hi 0 n —> oo T”T 

i—l 

(iii) For each i G { 1 , . . . , k} there exists a finite set of points S Cl 2 (which may 
be empty, but could contain up to k - 1 points) with the property that 



Unfan T 



in Wfoc(M?\S) as n — > oo. 



We refer to the map u^ : 5 2 S 2 as a ‘body’ map, and the maps uj 1 : 
S 2 — > S 2 as ‘bubble’ maps, or simply ‘bubbles.’ The points {x 1 , . . . ,x m } will be 
called ‘bubble points,’ and the A \ ‘bubble scales.’ Since each uj 1 is a nonconstant 
harmonic map between 2-spheres, the energy of each must be at least An. In 
anticipation of an application of Theorem (1.1) to the harmonic map flow, we now 
seek a generalisation of (1.3); we want to understand the behaviour of E for maps 
not necessarily C k - close to a harmonic map, but with small tension in L 2 . (Such 
maps need not be even W 1,2 -close to a harmonic map - see [11].) In [11, Theorem 
1.2] we proved such a generalisation, building on our earlier work in [10]. Here we 
state a weaker version with simpler hypotheses. 



Theorem 1.2. Suppose that we have a sequence of maps u n : S 2 — > S 2 satisfying 
the hypotheses of Theorem 1.1, and that we pass to a subsequence and find a limit 
Uoo, bubble points {x j } and bubble data uj 1 , >f n , a l n at each bubble point - as we 
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know we can from Theorem 1.1. Suppose that at each bubble point, only one bubble 
develops - i.e., that in the language of Theorem 1.1, k = 1 at each xT Suppose 
further that the energy density of the body map is nonzero at each bubble point - 
i.e., that |V^oo|(a: J ) ^ 0 for all j G { 1, . . . , ra}. Then there exist constants C > 0 
and k E N U {0} such that after passing to a subsequence, the energy is quantized 
according to 

\E(u n )-4irk\<C\\T(u n )\\i 2{s2) . (1.5) 

A further part of [11, Theorem 1.2] is a 4 repulsion’ estimate which controls 
some of the bubble scales by: 



Ki < ex P 



C||TK)|| 2 l2(52) J 



( 1 . 6 ) 



Remark 1.3. The full version of [11, Theorem 1.2] allows multiple bubbling at 
each point subject to a compatibility condition, and only requires the body map 
to have zero energy density at certain bubble points. It is at these bubble points 
that the bubble scales satisfy the repulsion estimate (1.6). 



One of the results of this paper is that both the quantization estimate (1.5) 
and the repulsion estimate (1.6) may fail if we allow arbitrary bubbles to occur at 
points where the energy density of Uoo is zero. See Corollary 2.2 below. Theorem 
1.1 gives us a basic understanding of the asymptotics of the harmonic map flow, 
as we now describe. By considering (1.2), given any sufficiently regular harmonic 
map flow u , we may find a sequence of times t n —> oo such that T(u(t n )) — * 0 in 
L 2 (S 2 ). We may then apply Theorem 1.1 to the sequence of maps (u(£ n )} to get 
bubbling convergence. (So far, this says nothing about the behaviour of the flow 
at times between t n and £ n +i0 It turns out to be possible (see [11] for details) 
to apply Theorem 1.2 to flows with appropriate bubbling, in order to obtain a 
quantization estimate for u(t) for all sufficiently large t - not just for the maps 
u(t n ). We state this result for smooth flows - note that the standard Struwe flow 
(as in [9]) is smooth for sufficiently large time, and we are only concerned here with 
the asymptotics of the flow. As before, the results in [11] are somewhat stronger. 

Theorem 1.4. Suppose that u : S 2 x [0, oo) — > S 2 is a smooth solution of (1.1). 
Then after picking times t n — > oo so that T(u(t n )) — > 0 in L 2 (S 2 ) and analysing 
the sequence {u(t n )} with Theorem 1.1, if only one bubble develops at each bubble 
point and the energy density of the body map Uoo is nonzero at each bubble point, 
then there exists C > 0 independent of t such that 

0 < E(u(t)) -E< C\\T(um\h { s 3), (1-7) 

for sufficiently large t, where 

E := lim E(u(t )) £ 47rZ. 

t—> oo 

Remark 1.5. Exactly as described earlier for smoothly converging flows, we can 
use an estimate such as (1.7) to deduce exponential decay of both [E(u(t)) — ~E\ 
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and f t °° ||T( , u(s))|| L 2 ( 52 )ds. This time, the consequence is the following corollary 
(cf. [11, Theorem 1.7]). 

Corollary 1.6. Under the hypotheses of Theorem 1.4, for any k G N and Q CC 
S 2 \{x 1 , . . . , x 171 } (where {x 1 , . . . , x 171 } are the bubble points) we have exponential 
decay of 

IK*) - Woo||l 2 (s 2 ) + IK*) - Woollen^)’ (1-8) 

as t — > oo. 

A priori, not even uniform convergence u{t) — » u ^ in L 2 (S 2 ) as t — > oo is 
clear. In fact, for more elaborate target manifolds than S' 2 , this may fail - see 
[10]. A further result of this paper is that if we allow a smooth flow to develop 
bubbles at points in the domain where the energy density of Uoq is zero, then the 
convergence of the flow need not be exponential - see Theorem 2.3 below. The 
combination of this result and [11] leads to apparently the first situation for the 
harmonic map flow in which the position of bubble points in the domain is shown 
to have a bearing on the properties of the bubbling. 



2. Details of results 

We first need the results of some energy and tension computations for a simple 
map between 2-spheres. 

Proposition 2.1. Given d G N and £ G (0,1], consider the map w £ : S 2 — * S 2 
defined in terms of stereographic complex coordinates by 

W e (z) =£ d (z d -Z~ d ). 

Then 

(a) code 2d < E(w e ) — 8nd < Code 2d , (2.1) 

for universal co,Co > 0; 

(b) Ci£ 2d_1 < ||T(w £ )|| L 2 (S 2) < Ci£ 2d_1 , (2.2) 

for Ci, Ci >0 dependent only on d. 

Given these calculations (see Section 3 for comments on the details) we may 

define a sequence of maps u n : S 2 — > S 2 by u n := Wi/ n to establish the following 

corollary. 

Corollary 2.2. For any d G N, there exist k G N and a sequence of smooth maps 
u n : S 2 — > S 2 with T(u n ) — » 0 in L 2 (S 2 ), such that 

2d 2d 

c\\T(u n )\\m 2) < \E(u n ) - Ank\ < C\\T{u n ) ||^ 2) , 

with c, C > 0 dependent only on d. In particular , the quantization estimate (1.5) 
fails for the sequence {u n }. 
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Moreover, if we apply Theorem 1.1 to the sequence {u n }, then the resulting 
bubble scales must satisfy 

c\\T(u n )\\ 

T 2 (S 2 ) — — C\\T{u n )\\ 

L 2 (S 2 ) ’ 

for some c, C > 0 independent of n. 

The second part of this corollary should be compared to the repulsion esti- 
mate in [11] mentioned in (1.6). The sequence u n converges to a constant body 
map with bubbles of degree d and —d. An alternative example, with slightly dif- 
ferent exponents, would arise by defining instead u n (z) = Wi/ n (nz). In this case 
we would get a degree d body map with a single degree —d bubble. Next we state 
precisely our result for the harmonic map flow. 

Theorem 2.3. Fix an integer d > 1 and define uo : S 2 — > S 2 by 

uq(z ) = z d — z~ d . 

Then there exist 77 = 77(d) > 0 and a unique smooth solution u : S 2 x [0, 00 ) — > S 2 
of the harmonic map heat equation (1.1) with 7/(0) = uq, and fort > 1 the solution 
has the following properties: 



(a) 


|| u(t) - U 00 || £ ,2 (S 2 ) > r)t 2(rf-i); 


(2.3) 


(b) 


\\T(u{t))\\ L 2 ( s 2 ) > 


(2.4) 


(c) 


J \\T{u(s))\\ms 2 )ds >rit <*-i; 


(2.5) 


(d) 


E{u(t)) — lim E(u(s)) > 77 


(2.6) 



Inequality (2.3) implies that u(t) does not converge exponentially to 
and should be compared to (1.8). (Note that for the flow in question, u(t) does 
nevertheless converge to Uqq as t — >• 00 .) Inequalities (2.5) and (2.6) imply in 
particular that their left-hand sides cannot decay exponentially to zero as t — > 00 , 
and should be compared to Remark 1.5. 

3. Elements of the calculations and proofs 

Although it is the complex coordinates which help us to pick natural maps for this 
problem, we are better off calculating in spherical polar coordinates (0, (/>) in order 
to later apply the comparison principle and exploit rotational symmetry. With 
our convention, 0 £ [0, n\ represents the ‘latitudinal’ angle from the vertical , and 
<fi measures longitude. We will be considering, for integral d > 1, maps of the form 

(M)^(/?(0),d0), (3.1) 

for appropriate (3 : [0, tt\ — > M. The energy of such a map is given by 

E = tt ( (P'(6)) 2 sin# + - 7 —z sin 2 (3(0) d0, (3.2) 

Jo L sin ^ J 
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and the tension may be written 



t = ^ 



where r(/3) : [0, tt] may be calculated 2 to be 

1 aun\ d 2 sin (3(9) cos 0(0) 
——0 ( 9 ) 



t(0) = 0"(9) + ~^0'(O) - “ (3 . 3) 

tan 0 sin 0 

abused notation by using (3 both as a function and as a target coor- 



and we have 
dinate. 

Proof. (Proposition 2.1.) The map w £ of Proposition 2.1 is of the form (3.1) with 
(3 = a where 



0 



a(0) = a £ (0) := 2 tan 1 e d ^tan d ^ - tan d ^ J . (3.4) 

Plugging this into (3.2) and using the substitution t = tan|, we eventually find 
that 



E-S7rd = 32nd 2 e 2d 



I 



+2d-l 



(1 + e 2d (t d -t~ d ) 2 ) 2 



dt , 



jo ' c ° ) J 

which then gives (2.1) of the proposition since the integral can easily be checked 
to be of order 1 uniformly for e E (0,1]. From the definitions of a and r, a 

Vitfrkrwarrl tvnf Inner ralmlatinn t.pilcps ns to 



tv L/V VI V1UVI ^ U1111V1 11X1J iv/l O j . i A V/J 

straightforward but long calculation takes us to 

—4 e M d 2 sec 4 f (tan rf ~ 2 f - tan~ d ~ 2 f ) 
(1 + £ 2d ( tan d | — tan -d f ) 2 ) 2 



(3.5) 



Plugging this into the expression 

\\'^( W e)\\ 2 L' 2 (S‘ 2 ) = 27T / 

Jo 

leads us to (2.2) of the proposition, with the integral being dominated now by 
nrrlcir nf tr n 



|r| 2 sin 0d0 , 



values of 0 of the order of 5. 



(3.6) 

by 

□ 



Remark 3.1. Related to w £ is the map u £ : S 2 — ► S 2 given by 

U £ {Z) =£ d (z d + Z- d ). 

Calculations for u £ show (cf. Proposition 2.1) that 

-c 0 d£ 2d < E(u £ ) - 8tt d < -C 0 de 2d (3.7) 

for universal Co, Co > 0, and 

CIS 2 *- 1 < ||T( We )|| L2(s2) < C^- 1 (3.8) 

for ci , Ci > 0 dependent only on d. This turns out to be useful in the theory of 
biharmonic maps, as may be described in [2]. 



2 cf. the case d = 1 in, say, [1] 
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We now consider the significance of the maps w £ to the harmonic map heat 
flow. As we see in the following proposition, it is possible to use these maps, with 
e dependent on time, as barriers. 

Proposition 3.2. Let e : [0,oo) — ► (0, 1] solve the ODE e — — 1 6d£ 2d ~ x with initial 
condition 6r(0) = 1. Then 

r\ 

-a e(t) {e)-T{a e(t) ){6)> 0 (3.9) 

for 9 E [0, §], t > 0. 

By symmetry, the left-hand side of (3.9) must be weakly negative for 0 E 
[|,7r]. The proof of the proposition is a straightforward calculation using the in- 
equality r 2 (l + s 2d (r d — r~ d ) 2 ) > e 2 /2 valid for r > 0, e E (0, 1], d > 1, say. Note 
that neither this inequality nor the ODE of the proposition is sharp. 

Remark 3.3. If we now have a smooth solution of (1.1) of the form 

(0,<M) -► (3.10) 

where h : [0,7r] x [0,T] — > M satisfies h( 0, •) = — 7r, h(0 , •) = —h(7r — 0 , •) and 
h(0, 0) < a\ (6) for 6 E [0, then by virtue of Proposition 3.2, we may apply the 
comparison principle (see [1] for details in a very similar situation) to establish 
that 

h(0,t) < a £{t) (6) 

for 0 E [0, and t E [0,T], with an analogous (reversed) inequality for 0 E [f , tt] . 
Similarly, if h{0, 0) E [— 7r,7r] for 6 E [0,7r], then the comparison principle tells us 
that h(0,t) E [— 7r,7r] for 0 E [0,7r] and t E [0,T]. 

Proof. (Theorem 2.3.) The work of Struwe [9] gives the existence of a weak solution 
u with u( 0) = uo, which is smooth on some maximal time interval [0,T) with 
T E (0, oo], and indeed smooth throughout space-time except possibly at finitely 
many points. There is no other smooth solution on [0,T) with uq as initial map 
- see for example [9]. Since the initial map is of the form (3.1) with f3 = aq, 
and ai(0) = — aq(7r — 6 ), the uniqueness we have described guarantees that the 
flow will be of the form (3.10) for t E [0,T), with h satisfying h( 0, •) = — n and 
h(8 , •) = —h(7r—9, •) for 6 E [0, \ ]. By the finiteness of the singular set, a flow of this 
form must be smooth away from the poles 6 — 0 and 6 = n at time t = T if T < oo, 
and so by the bubbling description of singularities in the Struwe flow given in [9] , 
the only way that u could fail to be smooth would be by bubbling at the poles. 
Meanwhile, since ft(-,0) = aq, Remark 3.3 tells us that — n < h(0,t) < a £ ^(0) 
for 6 E [0, and t E [0,T), where e(t) > 0 is as in Proposition 3.2. In particular, 
this restricts the image of the flow near the pole 6 = 0 at finite time, ruling out 
the development of a bubble at finite time and forcing T = oo. Note that our 
knowledge of the harmonic maps between 2-spheres, and the bubbling theory of 
Struwe then forces the flow limit Uoq to be the constant map ($,</>) — > (0 ,d<f). 
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Although this argument holds for all integral d > 1, we are assuming d > 1 in the 
theorem, and therefore the solution of the ODE for e in Proposition 3.2 is 

e(t) — (1 + 32 d(d — 1 )t)~ 2 ^ d -^ . 

For t > 1, say, we have 

e(t) > rjt ~ 2 ( d ~ 1 ) , (3-11) 

for some 77 = 77 (d) > 0 which we may reduce during the proof. Meanwhile, for 
0 < 6 < e < 1, we may estimate 

<* £ ( 0 ) < 

for some small 77 > 0 independent of £, and therefore, since a £ ^(0) > h(0,t) for 
this range of 6 , and all t > 0 , we have | u(t) - Uoo\ > 77 > 0 (for a possibly smaller 
77 ) on the disc defined by 6 < e(t). Consequently, 

II^W - Uoo\\ L 2 (5 2 ) > ^(t), 

which once combined with (3.11) is (2.3) of Theorem 2.3. In order to tackle (2.4) 
let us return to the definition of r given in (3.3), and assume that (3(0) = —it and 
= 0. Multiplying this definition by sin 2 0 (3'(0) and integrating from 6 = 0 to 
6 = ^ , we obtain (after integrating by parts) 

J\m 2 er(3\e)d6= l -{I3'(^)) 2 . 

Recalling (3.2) and (3.6), Cauchy-Schwarz then tells us that 

T\\l*- 

When we apply this to flow u and its corresponding h, keeping in mind the rela- 
tionship of h(*,t) and a £ ^, we find that 

l|2'(«(<))llL a (S’») > V (^(f>*)) ^ v ( a 'e(t)(^)) >ve 2d , 

where the value of 77 = 77 (d) > 0 is allowed to decrease at each instance. This estab- 
lishes (2.4). The final two parts of Theorem 2.3 follow easily from (2.4). Straight- 
forward integration gives (2.5), whilst integrating the square of (2.4), keeping in 
mind ( 1 . 2 ), gives ( 2 . 6 ). □ 

Remark 3.4. Only in the case d > 1 does the solution e(t) of the ODE in Propo- 
sition 3.2 fail to decay exponentially as required for Theorem 2.3. However, at £ (t) 
serves as a natural barrier even for d = 1 , preventing finite time blow-up - cf. 
Chang-Ding [1]. 



sin Or (3 f (0)d0 < CE*\ 
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Heat Flow into Spheres for a Class of Energies 

Norbert Hunger biihler 



Abstract. Let M and N be compact smooth Riemannian manifolds without 
boundaries. Then, for a map u : M — > N, we consider a class of energies 
which includes the popular Dirichlet energy and the more general p-energy. 
Geometric or physical questions motivate to investigate the critical points of 
such an energy or the corresponding heat flow problem. In the case of the 
Dirichlet energy, the heat flow problem has been intensively studied and is 
well understood by now. However, it has turned out that the case of the p- 
energy (p / 2) is much more difficult in many respects. We give a survey 
of the known results for the p-harmonic flow and indicate how these results 
can be extended to a larger class of energy types by using Young measure 
techniques which have recently been developed for quasilinear problems. 



1. Introduction 

l. 1. The p-energy 

Let M and N be smooth Riemannian manifolds without boundary of dimension 
m and n respectively. M and N are equipped with Riemannian metrics 7 and g. 
The p-energy of a differentiable map u : M — > N is defined to be 

E p (u) := - f \Du\ p dfi. 

P J M 

Here, Du : T X M — > T U ^N is the differential of u to which we associate the 
p-energy density 

e p (u) -\Du\ p = -(trac e((Du)* Du)) p ^ 2 , 

where denotes the adjoint with respect to the underlying inner product on the 
respective tangent spaces T X M and T U ^N. And, of course, dfi is the measure on 
M associated with the metric on M. In local coordinates, the p-energy can be 
expressed by the following formula: 

E p (u) '■=- [ o u)d Q u l dpu J ) p/2 y/j dx . 

P J M 

Here, = yj \ det(7 a /3)|, and we adopt the usual summation convention, i.e., we 
automatically sum over Latin indices from 1 to n and over Greek indices from 1 to 

m. By Nash’s embedding theorem, we always may assume that N is isometrically 
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embedded in some Euclidean space R k . In that case, the p-energy density is given 

by 

-(7 0/3 d a u ■ dpu ) p/ 2 
P 

where u : M — > N C R k is viewed as a map into R k and is the Euclidean inner 
product in R k . 

We remark, that for p = 2 the p-energy coincides with the well-known Dirich- 
let energy. Moreover, as it is easy to see, the p-energy is conformally invariant iff 
m = p. This fact allows to switch between conformal charts. 



1.2. The first variation of the p-energy 



Given a class F of mappings from M to N defined by boundary conditions on 
dM (if there is a boundary) and possibly additional topological conditions, e.g., a 
homotopy class, one may try to minimize the p-energy within the class F. One 
may hope that the solution of this minimization problem is a map u in F with 
particular analytic and geometric properties. In physics, the p-energy for maps 
into the sphere or into a Grassmannian manifold may serve as a model for liquid 
crystals. Stationary points of the energy then correspond to physical equilibria of 
the liquid. 

Necessarily, in a minimum u, the first variation of the p-energy must vanish 
at u for all variations of u with compact support in the interior of M. Let us 
consider variations u t : M — > N of u = uo for \t\ small such that u t (x) is of 
class C 1 in the variables (x,£) and u t (x) = u(x) for all t and all x outside some 
compact subset K C M which is contained in the domain U of a coordinate system 
x 1 , . . . , x 171 on M with range ft and is mapped into the domain of a fixed coordinate 
system y 1 , . . . ,p m on iV by all u t . Now, the support of the initial vector field of 
the variation 



V(x) := 



G T u (x)N 

t = 0 



is contained in K. We also assume that V is of the class C 1 . It is no restriction 
to assume u\(x) — u l (x) T tV l (x). If u is C 2 , the usual computation yields that 
vanishing first variation is equivalent to 



y/y dx@ \ lJ dx a dx$ ) dx a J 

V 13 dx a dx@ ) dx a dx@ 



(1) 



in U for Z = 1, . . . , n. Here, we used the Christoffel symbols 



r ij — 2 9 (9ik,j 9ij,k H~ 9jk,i ) 
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of the metric g with respect to the coordinates chosen on N. Thus, the right-hand 
side can be interpreted as 

A(u)(Du , Du)(pe p (u)) 1 ~ 2 ^ p 

where A(u)(-, •) denotes the second fundamental form of N. The operator on the 
left-hand side of (1) is called p-Laplace operator related to the manifolds M and 
N and is denoted by A p , i.e., A 2 is simply the Laplace- Beltrami operator of the 
manifold M and does not depend on N. A 2 is a linear elliptic diagonal operator 
in divergence form. For p > 2 (0 < p < 2) the operator is degenerate (singular) at 
points Du = 0. 

The right-hand side of (1) is for p = 2 a quadratic form in the first derivatives 
of u with coefficients depending on u. These strong nonlinearities are caused by 
the non-Euclidean structure of the target manifold N and cannot be removed by 
special choices of coordinates on N unless N is locally isometric to Euclidean space 
R n . But even in this case the space of mappings from M to N does not possess a 
natural linear structure unless N itself is a linear space. In general, the right-hand 
side of (1) is of the order of the p- th power of the gradient of u. 

A C 2 -solution u : M — > N of (1) is called p-harmonic map. 

Using variations in the domain , i.e., variations of the form 

u t {x) = u(x + t((x)) 

where £ = (C 1 , . . . , C m ) with each G C^(B p (y)), we obtain from ~u t (x) | t _ Q — 0 
an equation whose classical solutions contain the set of the p-harmonic maps. 
However its weak solutions, which are called “stationary p-harmonic maps”, need 
not contain the set of weak solutions of (1). In particular, p-energy minimizing 
maps are both, weakly and stationary p-harmonic. Weakly p-harmonic maps admit 
in general far worse singularities than the energy minimizing maps. Many results 
are known about stationary p-harmonic maps. See, e.g., Duzaar and Fuchs [9], 
Fuchs [18] and [19], as well as Strzelecki [41] and [40]. 

If we think of N as being isometrically embedded in some R fc , we consider the 
following variations: Let S C be a tubular neighborhood of N and tin : S — > iV 
the (smooth) nearest-neighbor projection. Denote T P N C TpR k the tangent space 
to N at a point p G N. Let <p G Cq(M, R k ) satisfy 

V(x) £ T U { X )N 

for all x G M and <p having compact support in a single coordinate chart of M. 
Then ip induces a Covariation at u : M — » N C 



u t = ttn o (u + tip ) . 
Now, the first variation of the p-energy is given by 



dt 



E p (ut) 



= / 

t = 0 J h 



ct(3 



du du \ 2 



M \ dx a dx@ 

/ A p u • ip -y/7 dx 

J M 



S-i 



q/3 



du dp 
dx a dx@ 



yfydx = 
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where A p is the p-Laplace operator related to M and R fc , i.e., 



A p u = 



_i d_ 

y/j dX@ 



/_ ( ^ a0 du j du j 
^ V dx° dx? 



P - 1 

o - L 



dx a 



Thus, the first variation vanishes at a C 2 map u if and only if 

A p u _L T U N . (2) 

We can make (2) more explicit by introducing a local orthonormal frame v n + 1 , 
. . . ,Vk for ( T P N ) ± , the orthogonal complement of T p N in R k . Then, by (2) there 
exist scalar functions A n+1 , . . . , X k such that 

k 

-A p u = A l (yi o u ) . (3) 

l=n -\- 1 



Multiplying (3) by i/j o u (i fixed) and using the fact that - — • ui (u) = 0 for all 

UX a x ' 

a, we obtain 



A* = 



r 



aP <^dw>_\ 2 1 _, a0 dip duj (n) 



dx a dx@ J ^ dx a dx@ 

The Euler-Lagrange equations of the p-energy can also be formulated in the fol- 
lowing coordinate free form: 



( 4 ) 



e{u)r{u) + (1 



2 

P 



)Dugr&d 1 e(u) — 0 . 



( 5 ) 



Here, r(u) denotes the tension field of u : r(u) = trace 7 Vfiw, and V denotes the 
pull-back covariant derivative in the bundle T*M<g>u~ 1 TN. For p = 2 the equation 
reduces simply to 

t(u) = 0 . 

As an application, we see, that the identity mapping id m : M — > M of a Riemann- 
ian manifold is p-harmonic: Since T>(idM) has constant coefficients with respect to 
all coordinate systems on M, VD(idM) = 0 and the first term in (5) vanishes. On 
the other-hand, the p-energy density e(idAf) is constant and hence the gradient in 
the second term disappears, too. We remark that in general idM is not p-energy 
minimizing within its homotopy class (see Eells-Lemaire [11]). 

The natural space to work in when considering questions concerning the p- 
energy is the nonlinear Sobolev space 



W 1,P (M, N) := {/ G W hp (M,R k ); f{x) e N for //-almost all x G M}. 



Notice, that if both M and N are compact, then different embeddings of N in R k 
give rise to homeomorphic spaces W 1 ^ P (M,N). For an intrinsic definition of the 
space W 1,P (M, N) see Federer [17]. If N is isometrically embedded in R fc , then 
the formulas defining the p-energy for C 1 functions make sense also for functions 
in W 1,P (M, N). If a function u G W 1,P (M, N) is a weak solution of the Euler- 
Lagrange equations we discussed previously, then it is called a weakly p-harmonic 
map. 
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The space FT 1,P (M, N) defined as the closure of the class of smooth func- 
tions from M to N in the VF 1,p -norm is contained in W 1 ^ P (M, N) but does not 
coincide with the latter space in general (this fact gives rise to the so called “gap 
phenomenon” of Hardt-Lin [21]). This important observation was first made by 
Schoen and Uhlenbeck: see Eells and Lemaire [11] as a main reference. However, we 
have iL 1,p (M, N) = FF 1,P (M, N) if dim(M) = p (see Schoen and Uhlenbeck [39], 
Bethuel [2] or Bethuel and Zheng [3] ) . The gap phenomenon limits the possibilities 
to approximate the p- harmonic flow, which we will discuss next, by solutions with 
smooth initial data. 

1.3. The p-harmonic flow 

A basic existence problem for p-harmonic maps is the homotopy problem: 

Given a map uo : M — > N is there a p-harmonic map u homotopic to uo ? 

This question was first answered in the case p = 2 of harmonic maps. The 
answer is affirmative if the sectional curvature K N of N is non-positive (see 
Eells-Sampson [12]), or - in case of a two-dimensional surface M - if the sec- 
ond fundamental group of N is trivial: ^(N) = 0 (see Lemaire [30] and Sacks- 
Uhlenbeck [38]). Eells and Wood destroyed the hope for a more general theorem 
by the following counterexample in [13] by showing that, if u : T 2 — > S 2 is 2- 
harmonic, then deg u ^ ±1. Another counterexample has been given by Lemaire 
in [30]. Thus, in general the attempt to solve the homotopy problem by minimiz- 
ing E within a given homotopy class will fail: Homotopy classes are not weakly 
closed in W 1 ' P (M, N) generally. Thus, it may be difficult or impossible to solve 
the homotopy problem for p-harmonic maps by direct variational methods. As the 
key idea to get around this difficulty, Eells and Sampson proposed in [12] to study 
the heat flow related to the 2-energy: 

ut — A 2 ^ = A(ix)(Vu, Vu)m on M x [0, oo[ (6) 

with initial and boundary data 

u — uo at t = 0 and on dM x [0, oo[ 

for maps u : M x [0, oo[— > N C Here, A(u) : T U N x T U N — ► ( T U N )-*- is the 
second fundamental form of N. The idea behind this strategy is of course that a 
continuous deformation u(-,t) of uo will remain within the given homotopy class. 
Since (6) may be interpreted as the L 2 -gradient flow for the 2-energy, one may 
hope that the solution u(-,t) for t — » oo will come to a rest at some critical point 
of E2 that is a harmonic map. For target manifold AT, satisfying the geometric 
restrictions mentioned above, this program has been applied with success for p = 2. 

Another approach to the homotopy problem for p-harmonic maps has been 
given by Duzaar and Fuchs in [10]: they extended the Eells-Sampson result to 
the case p G [ 2 , 00) by using an asymptotic analysis of the not degenerate energy 
/ m ( £ + \ D u\ 2 ) p/2 dfi, S> 0 . 

Corresponding to the harmonic flow, the p-harmonic flow is described by 
u t — A p u = (pe(u)) 1 ~p A(u)(Du, Du) on M x [0,oo[ 



(7) 
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with initial data 

u = uo at t = 0 

for maps u : M x [0, oo[— ► N C R k . 

We now want to briefly describe the known results concerning the p-harmonic 
flow. We restrict this presentation to the case p ^ 2. 

The p-harmonic flow was first considered in [6] and [22]. There, the following 
result has been proved for the p-harmonic flow into spheres, i.e., for 



dtu — A p u — pe p (u)u 


on M x [0, 00 [ 


(8) 


u(-,0) = u 0 


on M 


(9) 


\u\ = 1 


pi- a.e. on M x [0, 00 [. 


(10) 



Theorem 1. For initial data uq E W 1,p (M, S n ), p >2, there exists a global weak 
solution u to the equation (8) -(10). This solution is weakly continuous in t > 0 
with values in W 1,P (M), i.e., for any test function g E C°°(M), hi(t) = J M u-gdp, 
and h 2 (t) = f M Du • Dgdp, are in C 0, 2 (R + ). Furthermore for almost every t > 0 
this solution satisfies the energy inequality 

[ \\9tu\\l HM) dt + E p (u(t)) < E p {u 0 ) . (11) 

J 0 

The technique was the same as in the corresponding case of the harmonic 
flow which was solved by Chen in [5]: The idea is a penalization technique to 
approximate the p-harmonic flow, which will be explained and used in Section 2 
for a more general class of energies. The monotonicity of the p- Laplace operator 
allows to solve the approximating equations by Minty’s trick alone. Nonetheless, 
to pass to the limit in the penalized equations some stronger compactness results 
on the p-Laplacien are needed. In the original work, this is based on the uniform 
monotonicity of the p-Laplace operator: 

Theorem 2. For k = 1,2,..., let fk : M x [0, T] — > R z satisfy the equation 
d t fk ~ \fk = gk, on M x [0, T] 
in the sense of distributions. Assume that 

1. {fk}keN is bounded in L°°(0,T; W 1,P (M, R 1 )), 

2. {dt fk}ke n is bounded in L 2 (0, T; L 2 (M, R z )), and that 

3. {gk}ke n is bounded in L l {0 , T; L 2 (M, M z )). 

Then, {fk}keN is precompact in L 9 (0,T; W 1 ^ q (M, R 1 )) for each 1 < q < p. 

In Section 2 we offer a new approach to this compactness problem which 
involves some Young measure techniques, and which does not use strict or even 
uniform monotonicity properties of the operator. This allows in particular to obtain 
results for the flow of energies which are merely convex. 

Theorem 1 was later proved for 1 < p < 2 by Liu in [32], [33] and by 
Misawa in [37]. The p-harmonic flow from a unit ball in M m into S' 1 C R 2 was 
also considered by Courilleau and Demengel in [7]. There, also non-uniqueness of 
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the p-harmonic flow with weakly p-harmonic, but not stationary harmonic initial 
data was proven, a result, however, which is already contained in [44] and even 
before that in [25]. 

The p-harmonic flow equation has two nonlinear terms: The p-Laplace term, 
and the right-hand side involving the second fundamental form. The previously 
described compactness techniques allow to pass to the limit in the p- Laplace term 
only. The big problem, which is unsolved to the present day, is to pass to the 
limit in the nonlinear term on the right-hand side of the equation. The sphere as 
target manifold has enough symmetries to rewrite the equations in such a way, 
that the right-hand side simply vanishes (compare Section 2). A similar trick can 
be played to prove existence of the p-harmonic flow if the target manifold A" is a 
homogeneous space with a left invariant metric. 

Theorem 3 ([23]). For 2 < p there exists a global weak solution of the p-harmonic 
flow between Riemannian manifolds M and N for arbitrary initial data having 
finite p- energy in the case when the target N is a homogeneous space with a left 
invariant metric. The solution u: M x [0, oo[— » N satisfies the energy inequality 

f \d t u\ 2 dtdn + - [ \Du(T)\ p d/j,<- [ \Dum p d/j, ( 12 ) 

M P JM P JM 

for almost all T > 0. 

In this case of a homogeneous space as target, the equation of the p-harmonic 
flow can be reformulated such that the right-hand side has a divergence structure 
(see [24] and [23]). In [23], a technique different from the penalization technique 
was used to approximate the equation, namely a time discretization. The reason 
is, that the penalization does not preserve the special structure of N. 

The analogous result for 1 < p < 2 was proven in [31]. There, the case of 
a modified p-energy is treated, namely that of the p-energy ‘with potential’, i.e., 
E{u) == E p — f M H(u)dpi for a Lipshitz function H. 

Another particular case which has been solved is the conformal case p = m = 
dim M: 

Theorem 4 ([26]). For given initial value uo € W 1,m (M, N) there exists a weak 
solution u : M x [ 0, oo[— » N of the m-harmonic flow 

u t - A mu _L T U N 
u(-,0) = u 0 . 

u satisfies the energy inequality (12) and is in W 1,rn (M ) weakly continuous in 
time. There exists a set E = U^ =1 E^ x {Tk}, E/e C M, 0 < Tk < oo, such that 
on every open set Cl C M x [0, oo[ with dist(fi, (M x {0}) U E) = p > 0 there 
holds ||VT/||co,/3(n) < C for some constants C (depending on m, E m (uo), M, N 
and i a) and (3 e]0, 1[ (depending on m, M and N). There exists £ i > 0 such that 
the number K of singular times is a priori bounded in terms of the initial energy, 
K < £f 1 E rn (uo), and the singular points (x,Tk) are characterized by the condition 
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limsup^ Tfc E^uft), Br(x)) > £i for any R > 0. At every singular time T k the 
decrease of the m-energy is at least £i: 

E m (u(T k )) < lim inf E m (u(t)) - ei. 

t/STk 

In order to get rid of the nonlinear side condition, that the solution has to 
take values in the sphere, Hamilton’s geometric technique from [20] was used, 
namely that of a total geodesic embedding of the target N in some R k . I.e., the 
Euclidean metric in R k is deformed to a new metric h such that 

1. N is still embedded isometrically, 

2 . the metric h equals the Euclidean metric outside a large ball, 

3. there exists an involutive isometry t : T — > T on a tubular neighborhood T 
of N corresponding to multiplication by -1 in the orthonormal fibers of N 
and having precisely N for its fixed point set. 

A h- geodesic curve 7 connecting x, y G N (x, y close enough) will always be 
contained in N. This follows from the (local) uniqueness of geodesics and the fact 
that with 7 the curve t o 7 is another geodesic joining x and y. 

The idea is that, if the initial data are smooth, then a short time solution u 
of the p-harmonic flow is unique which forces it to stay in N (since otherwise, u 
and lou were two different solutions) . 

The crucial point in the conformal case is an energy concentration lemma: 

Theorem 5 ([26]). If dim(M) = p and u G C 2 (M x [0,T[;7V) is a solution of 
the m-harmonic flow then there exist constants c, e 0 > 0 which only depend on 
the geometry of the manifolds M and N, and there exists a time T 0 g]0,T[ which 
depends in addition on uo, with the following properties: If the initial local energy 
satisfies 

sup E rn (u(0),B 2R (x)) < £ 0 

x£M 

then it follows 

E m (u(t), Br(x)) < E m (u(0), B 2 r(x)) +cE 1 ^ 1 ‘E (13) 

for all (x,t) G M x [0, To]. Here E 0 denotes the initial energy, and E m (u(t),B) is 
the m-energy of u(t) in a region B C M. 

This lemma allows to control energy concentration of the flow for short time 
which leads to short time existence and (local) smoothness. At times, where the 
energy concentrates, the flow develops a singularity and can be restarted. It is 
not known whether for m /2 the m-harmonic flow develops singularities in finite 
time, in contrast to the harmonic flow, where this is confirmed by the example of 
Chang, Ding and Ye in [4] . 

Technically, the proof of Theorem 4 is carried out by two approximation 
steps: In a first step, only smooth initial data are considered. At the same time, 
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the p-energy is regularized by E p , £ = f M {e + \Du\ 2 ) p ^ 2 dpL. In a second step, arbi- 
trary initial data are approximated by smooth functions (observe, that by Bethuel- 
Zheng [3] C°°(M, N) is dense in W 1,P (M, N)). Misawa improved Theorem 4 in [34] 
by showing, that, as it is the case for the harmonic flow, the energy concentration 
set E is discrete not only in time, but also in space. 

The p-harmonic flow in the case of non-positive sectional curvature of the 
target manifold was treated by Fardoun and Regbaoui in [14] and [16]: 

Theorem 6. If Riem^ < 0 and uo G C 2,a (M, N), then forp > 1: 

1. There exists a unique global weak solution u of the p-harmonic flow from M 
to N such that d t u G L 2 ([0, oo[xM) and u,Du G C^([0, oo[xM) for some 
[3 g]0, 1 [. Moreover, the solution satisfies the energy inequality (12) for all 

t > 0. 

2. There exists a sequence tk — > oo such that u(tk , •) converges in C 1 ^ (M,N), 
for all /3 f < (3, to a weakly p-harmonic map Uoo G C 1,/3 (M, N). 

In particular, this results settles the homotopy question in this case. The line 
of proof goes along regularization of the p-energy, and a Bochner formula for the 
energy density (see [14] for details). 

The same result for p > 2 was shown earlier by Misawa in [35] under the 
assumption that the image of u$ is contained in a geodesic ball, and in [36] without 
that latter assumption. 

The p-harmonic flow for the case of small initial data has been treated by 
Fardoun and Regbaoui in [15]: 

Theorem 7. For any p > 1 and K > 0 there exists an e o > 0 depending on 
K,M,N and p such that if uo G C 2,a (M, N), 0 < a < 1, with E p (uo ) < £o and 
\\Du 0 \\ L°° (M) < K, the following is true: 

1. There exists a unique global weak solution u of the p-harmonic flow from M 
to N such that d t u G L 2 ([0, oo[xM) and u,Du G C^([0, oo[xM) for some 
(3 G]0,1[. Moreover, the solution satisfies the energy inequality (12) for all 
t> 0. 

2. There exists a sequence tk — > oo such that u(tk , •) converges in C 1, P ( M,N ), 
for all (3 r < (3, to a weakly p-harmonic map G C 1 ^ (M, N ) . Moreover 
there exists £o > 0 depending on K , M, N and p such that if in addition 
E p (uq) < £o, then Uoq is a constant map. 

This result of Fardoun and Regbaoui concludes our survey on the p-harmonic 

flow. 



2. Heat flow into spheres for a class of energies 

In this section, our target manifold N is the unit sphere S n C M n+1 . In order to 
keep the formulas short, we restrict this presentation to the case M = T = M m /Z m 
of a flat torus as domain manifold. The modifications necessary for a general 
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compact smooth Riemannian manifold without boundary are purely technical and 
straightforward. 

The class of energy functionals we want to investigate here is of the form 

E(u) = L p(\Du\)dx 

for functions u : T — » S n C M n+1 . Here p : M> 0 — > M is supposed to be continuously 
differentiable, convex, and satisfies the coercivity and growth conditions 

< p(0 for all f > 0 (14) 

0 < p'( 0 < C 2 e~ l for all £ > 0 (15) 

for two positive constants C\ and C 2 , and for a given p > 2. The energy flow 
related to this energy is described by the following equation 

d t u -div {p\\Du\)-~^ J_ T U S. 

We use the shorthand notation 

a(Du) := p \\Du\)^- 
I Du\ 

and hence the previous equation can be rewritten in the form 

d t u — div cr(Du) = Xu 

for a function A (x,t). Observe, that by multiplying this equation by u , we find 

A = cr(Du) : Du. 

Here, A : B := trace(^4H T ) denotes the usual inner product of two matrices of the 
same type. So, the final form of the heat equation we want to solve, is 

d t u — di va(Du) = ucr(Du) : Du on T x M >0 (16) 

u(-,0)=u o on T. (17) 

We suppose, that the initial data u 0 have finite energy, i.e., u 0 £ ll /l p (T, S n ). 

2.1. Approximation by penalization 

We construct a family of approximating problems by the classical penalization 
method. To this end, we consider the energies 

E k (u) := E{u) + k f x(M )dx 
Jt 

for functions u :T -> R n+1 , and for k £ N. The smooth function x : ®>o ^ M is 
chosen in such a way, that 0<x<l,x=l outside the interval [|, §], \{x) = 0 if 
and only if x — 1 , and that changes sign only once (see Figure 1). The idea is, 
that we abandon the restriction that u takes values in the target S n , but that (in 
terms of energy) with increasing k, it becomes more and more favorable for the 
function u to take values close to the sphere S n . 
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The heat flow equation for the penalized energy is 

d t u — di va(Du) + kx\\u\)^-r = 0. (18) 

M 

The goal is now to prove existence of a global weak solution for this penalized 
energy flow with the given initial data (IT). 

Lemma 8. Let k G N and uo G W l,p (T, S n ). Then, (17)-(18) possesses a weak 
solution Uk G L°°(0, oo; W l,p (T, R n+1 )) which is subject to the energy inequality 

[ \\dtUk\\ 2 L 2 (T) dt + Ek(uk(r)) < E(u 0 ) (19) 

Jo 

for almost all r > 0. 

Proof. W 1,P (T , M n+1 ) is separable and possesses therefore a Galerkin base {wi}^, 
Wi G W 1,P (T, E n+1 ) in the sense, that for every w G W 1,P (T, R n+1 ) there exist 
coefficients such that 



3 

w in W 1,P (T) as j — ► oo. 
i= 1 

It is convenient to choose the Wi smooth and L 2 -orthonormal. For fixed j G N, we 
make the ansatz 



3 

u^\t) Y^\t)wj 



2—1 

for a solution of the approximating system 
wi+cr(Du^) : Dwi + kx'(\u^\) 
with initial values 




\u(L\ 




= 0 



for l = l,...,j (20) 



i/^(0) = — > uo in W 1,P (T) for j — ► oo. (21) 

This is a system of j ordinary differential equations for the coefficients c^\t), 
l = 1, . . . , j. It is easy to verify that this system satisfies the hypotheses of the 
existence theorem of Picard- Lindelof and that therefore a solution exists on a time 
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interval ]0,r(j)[. Since, by construction, is an element of span(rci, . . . , Wj), 
equation (20) holds with in place of wi, and hence, after integration over a 
time interval ]0,r[, r < r(j), one gets 

i||^^^(r)||| 2 ( T ) + f f a(Du^) : Du^dxdt 

2 Jo Jt 

/»T /* (o') "1 

+ Jo j T kX '^ u(3) ^^W\ uU)dxdt = 2^° )( ' 0 ^^ 2 ( t ) ( 22 ) 

The first term in (22) is equal to l|c^(r)|jj^ . By (15), the second term is non- 
negative. Since \x'(t)\t < Imax^l for t > 0, the third term, in absolute value, 
is bounded by Cr for a constant C. This shows, that the functions c^\t) are 

bounded on [0,r(j)[, and hence, by (20), the same is true for Therefore, 

the existence interval is both open and closed and hence we obtain a global solution 
vfj) of (20). 

Similarly, with d t u^ in place of wi in (20), we obtain 

f \\d t u^\\ 2 L 2 ^dt + Ek(u^\r)) = Ek(u^) < C for all j and r > 0. (23) 

Jo 

Observe, that 

Ek(u ) > E(u) = f p(\Du\)dx > Ci f \Du\ p dx 

Jt Jt 

by the coercivity assumption (14) on p. Therefore, by (23), the sequence {u^}j 
is bounded in L°°(0, oo; W 1,P (T)) On the other hand, still by (23), the sequence 
{d t u^}j is bounded in L 2 (Tx]0, oo[). Aubin’s lemma therefore implies that there 
exists a subsequence {u^ ) }j> and a function u such that 

vJi } — > u in L r (Tx]0,r[) for all r e]0,oo[ and all r E [p, — — (24) 

m — p 

Notice that to have one subsequence which works for all such r and r, one can 
apply the usual diagonal sequence technique. By passing to a further subsequence 
if necessary, we can also assume that 



) — > u point- wise almost everywhere on Tx]0, oo[ 


(25) 


and 




U U') u weakly* in L°°(0, oo; W hp (T, R n+1 )) 


(26) 


and finally, that 




d t u^ d t u weakly in L 2 (Tx]0, oo[). 


(27) 


Notice also, that by the growth condition (15) we imposed on p, 
{cr (Du^)}j is bounded in L°°(0, oo; L p (T)). Therefore, 


the sequence 


a(Du ^ ^) w weakly* in L°°(0, oo; L p (T)) 


(28) 


for a further subsequence. 






Heat Flow into Spheres for a Class of Energies 



57 



For a test function </)(x,t) = £( t)rj(x ), with £ € Co°(]0, oc[) and 77 G 
span(ici, . . . , ic/), we infer from (20) that (as soon as J' > /) 



H dM j ' : 

Jo Jt 



4> + 



n a{Du {j,) ) : D4> 

I 



+ 



n /./x ii 
| as / — > oo 



[°°[ dM + r[w:Dcf> + =0. 

Jo Jt Jo Jt Jo Jt \ u \ 

(29) 

To pass to the limit in the first and the second term of (29), we used (27) and (28). 
For the third term, we used (25) and Lebesgue’s dominated convergence theorem. 
By density of the linear span of test functions we used, the resulting equation in 
(29) follows also for arbitrary </> G Co°(Tx]0, oo[). 

In order to identify w we now want to apply Minty’s lemma for monotone op- 
erators. (see, e.g., [42, “decisive monotonicity trick”]). To do this, we first observe, 
that by (22), for f — * oo, and by using (24), (25) and (27) 



f ( — <j(Du d )) : Du^ ^ = f f d t u^ ^ u ^ ^ -f f f { \u^ ^|) . r/ J 

Jo Jt Jo Jt Jo Jt \ uU } \ 

-> [ [ d t uu+ [ [ kx'(\u\)-Tu. (30) 

Jo Jt Jo Jt \ u \ 

On the other hand, by (29), we have 

[ [ -w : Du = f [ ud t u+ f ( kx'{\u\)^-u (31) 

Jo Jt Jo Jt Jo Jt \ u \ 

We now apply Minty’s lemma for the monotone operator 

A : L p (0,t; W 1,p (T)) -► z/(0, r; W~ l ' p ' (T)) 

u ■ — > ^ f J a (D u ) '■ ■ 

Now, as desired, from (26), (28), and the fact that the right-hand sides of (30) 
and (31) agree, it follows from Minty’s lemma that w = cr(Du). Hence, u is a weak 
solution of (18). 

The energy inequality (19) for u follows from (23) as j — > oo: We have 



,U') 



lim inf 

j-^oo 



[ \\dtu^ ||i 2 ( T ) > f \\dtu\\ 2 L 2( T j 
Jo Jo 



by (27), 

lim inf Ek{u^\r)) > Ek(u(r)) for almost all r 

oo 

by (26), the convexity of p, and by (25). And finally, observe that 

lim E k (ul 3) ) - E k (u 0 ) 

by (21), the Vitali convergence theorem and by (25). □ 
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2.2. Passage to the limit I 

We now want to pass to the limit in (18). Observe first, that by the energy in- 
equality (19), the sequence {uk} we constructed is bounded in L°°{ 0, oo; W 1 ,P (T)) 
and {d t Uk} is bounded in L 2 ( 0, oo; L 2 {T)). By passing to a subsequence which we 
still denote by Uk, we may assume that 

Uk u weakly* in L°°(0, oo; W 1,P (T)) (32) 



and 

d t Uk — ^ dtu weakly in L 2 (Tx]0, oo[). (33) 

Moreover, by Aubin’s Lemma, we can extract a further subsequence such that 

Uk — u in L r (Tx]0,r[) for all r e] 0, oo[ and all r E [p, — — (34) 

m — p 

and hence, for yet a further subsequence 



Uk — > u point- wise almost everywhere on Tx]0, oo[. (35) 

Notice, that from the energy inequality (19) it follows, that \uk\ — ► 1 almost 
everywhere. Combining this with (35), we conclude, that \u\ = 1 almost everywhere 
on T x [0, oo[, i.e., u takes values in the sphere, as desired. This is, nonetheless, 
not sufficient to already pass to the limit in (18). 

The following lemma is a maximum principle for the weak solutions Uk of (18). 



Lemma 9. Let Uk be a weak solution of (18) on [0,r[. Then, \uk(t)\ < 1 for all 
t E [0, r[, if\u k (0)\ < 1. 



Proof. To ease notation, we write u in place of Uk • We use 

u 

R(u) : — — u -j 7 minfl, ltd} 

\u\ 

as a test function in (18) and find for the three single terms: 

[ R{u)d t u = d t [ £(M) 

Jt Jt 



for 



Therefore 



0 



«*) = i 



i(x-l) 2 



if 0 < x < 1 
if x > 1. 



[ R(u)d t u = d t [ (\u\ - l) 2 

Jt J{\u\>i} 

In the second term, we get 

/ k X '(\u\)^-R(u) > 0 
Jt M 

since x'(\ u \) ^ 9 f° r 1^1 ^ 1 an d R(u) = 0 for \u\ < 1. For the third term, we have 

| u Du \ 2 

jyisu) . ; -r / fjy\±su\)- 



f a {Du) : DR(u ) = f a{Du ) : Du{ 1 - -*-) + f p f (\Du\) 

Jt J{\u\>i} \ u \ J\\u\>i} 



\Du\\u \ 3 
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Observe, that the last term in the previous expression is non- negative. Hence, 
summarizing, we obtain 



ld_ 
2 dt 




1)2 + l, 



1 



a {Du) : Du{ 1 — — ) < 0. 

i\u\>i} M 



Since the second term is non-negative, the derivative in the first term is non- 
positive, and therefore the lemma follows. □ 



By (31), we have for arbitrary r > 0 

/ / d t u k u k + / / a(Duk) ■ Du k + k / / x'(Wk\)r^-;U k = 0. 

Jo Jt Jo Jt Jo Jt \ u k\ 

The energy inequality (23) and the growth condition (15) imply, that the first two 
terms are, for fixed r, uniformly bounded in k. In view of Lemma 9, and the fact, 
that vanishes on [0, \] we obtain 




kM) W\ 



for a constant C which does not depend on k. 

We conclude, that the right-hand side in 

d t u k -di va(Duk) = -kx'(\u k \) t^t 

\Uk\ 



(36) 



is a bounded sequence in L X {T x [0,r]). This fact suggests now to consider the 
Young measure generated by the sequence Uk 0 Duk . Observe, that Uk converges 
almost everywhere, an therefore, if Duk generates the Young measure v on the 
space M( n+1 ) xm of (n + 1) x m matrices, then Uk 0 Duk generates the Young 
measure 5 U 0 v (see, e.g., [27]). We use the next section to collect some relevant 
facts about the Young measure v. 



2.3. The Young measure generated by Du & 

Lemma 10. Let v be the Young measure generated by the the gradients Duk of the 
solutions of (18) constructed in Lemma 8. Then 

(i) V( x ,t) i s a probability Young measure for almost every (x,t) E T x [0,r]. 

(ii) V{ x ,t) satisfies Du{x,t) — {u^ x t ),id) for almost every (x,t) E T x [0,r]. 



Proof (i) follows directly from the fact, that Duk is bounded in L l {T x [0, r]) (see, 
e.g., [1] or [27]). 

(ii): Duk converges weakly to Du in L P {T x [0, r]). On the other hand, 
since Duk is an equiintegrable sequence, we have that Duk converges weakly to 
in L l (T x [0,r]). Therefore, Du(x,t) = (i/ (a . )t) ,id). □ 



The final ingredient we need is a “div-curl inequality” : 






(er( A) — <j{Du)) : (A — Du)dv^ x t ^ (A) dxdt < 0 . 



(37) 
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This inequality appears first in [8] for the Young measure v generated by sequence 
of the gradients of solutions of a stationary non-linear elliptic system of p-Laplace 
type with right-hand side bounded in L l . The class of operators considered there 
includes in particular our case cr(Du). Due to the fact that, as we have seen, 
the right-hand side in (36) is bounded in L 1 , the proof in [8] carries over to the 
parabolic case. 

2.4. Passage to the limit II 

The aim is now, to pass to the limit in the penalized equation (18). The problem 
in the operator term is, that p need not be strictly convex, and hence a is not 
necessarily strictly monotone, or even uniformly monotone like in the case of the 
p-Laplace operator. Nonetheless, if p is not strictly convex in a region, it must be 
affine there, and it should therefore nevertheless be possible to pass to the limit 
there. It turns out that Young measures provide the adequate tool to go through 
with this program. 

We start with some notation. 

We consider the limit u of Uk as in (32)-(35). Let (x,t) be a fixed point in 
T x [0, r]. Then, let 

L := {(X,p(\Du(x,t)\) +a{Du(x,t)) : (A - Du(x,t))) \ A eM (n+1)xra } 

c M (n+l)xm x R 

denote the supporting hyperplane to the graph of p( |A|) in ( Du(x , t),p(\Du(x, t)|)) 
(see Figure 2). Moreover, let 

L 



A G M( n+1 ) xm 



Figure 2. Graph of p 

K {x , t) = {A € | p(|A|) = p(\Du(x,t)\)+a(Du(x,t)) : (A - Du(x,t))} 

denote the set, where p agrees with the supporting hyperplane L (see Figure 2). 

Lemma 11. For almost all (x,t) G T x [0 , r], the Young measure V{x,t) is supported 
in the set i.e spt iF(^x,t)) 
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Proof. By (37) the integrand (cr( A) — cr(Du(x,t))) : (A — Du(x,t)) (which is non- 
negative everywhere because of the monotonicity of cr) must vanish as a function 
of A on the support of the measure V( x ,t)- This is true with the possible exception 
of a set N of (x,£) G T x [0,r] of measure zero. Let us fix A G spt(^( x ^)) for 
(x, t) £ N, then 

(1 — a)(a(Du(x , t)) — <r(A)) : ( Du(x , t) — A) = 0 for all a G [0, 1]. (38) 

On the other hand, again by the monotonicity of a, we have 

0 < (1 — a)(cr(Du(x,t) + a(\ — Du(x,t))) — cr(X)) : (Du(x,t) — A) for all a G [0, 1]. 

(39) 

Subtracting (38) from (39) we obtain 

0 < (1 — a)(a(Du(x , t) + a(X — Du(x , £))) — a(Du(x , t))) : ( Du(x , t) — A) 

for all a G [0,1]. (40) 

But by monotonicity, also the opposite inequality holds true: 

0 > (1 — a)(a(Du(x , t) + a(X — Du(x , t))) — a(Du(x , t))) : (Du(x, t) — A) 

for alia G [0,1]. (41) 

Therefore, still for the same A G spt(i/( x?t )), by (40) and (41), we have 

0 = (1 — a)(cr(Du(x , t) + <a(A — Du(x, t))) — a(Du(x , £))) : ( Du(x , £) — A) 

for all a G [0, 1]. (42) 

Using (42), we conclude 

p(|A|) = p(\Du(x,t)\) + f cr(Du(x,t) + a(A — Du(x,t)) : (A — Du(x,t))da = 

Jo 

= p(\Du(x,t)\) + cr(Du(x,t)) : (A — Du(x,t)) 

and hence, A G K( x ,t)> as claimed. □ 

Lemma 12. (i) For each X G K( x ,t) there holds a(X) — <r(Du(x,t)). 

(ii) For almost all (x, t) G T x [0, r] following is true: 

[ a{\)dv {x ^ t) {\) = a(Du(x,t)). (43) 

J _A/f( n + !) x m 

Proof (i) follows the fact that p is continuously differentiable. 

(ii): As we have seen in Lemma 11, the support of the measure V( x ,t) ls 
contained in the set K( x ^y Therefore, in (43), we only need to integrate over 
K [Xit) C But there, as stated in (i), <r(A) = cr(Du(x,t)). Hence, the 

claim follows from the fact that, by Lemma 10(i), V( x ,t) ^ as mass 1. □ 

Actually, from (ii) in the previous lemma and the fundamental theorem on 
Young measures (see, e.g., [1]) it now already follows, that a(Duk) — ^ a (Du) in 
L l (T x [0 ,t]. However, as we will see in the next lemma, this convergence is even 
strong. 
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Lemma 13. Let Uk be the sequence solutions of (18) constructed in Lemma 8, u 
the limit of this sequence in (32), and 

9k(x,t) := \a(Duk(x,t )) - cr(Du(x, t))\ 

Then , gk — > 0 in L l (T x [0 , t]) for all r > 0. 

Proof Observe, that, since {Duk} is bounded in L°°(0, oo; L P (T)), by (15) cr(Duk) 
is bounded in L°°(0,r;L p (T)). Therefore, gk is an equiintegrable sequence on 
T x [0, t]. Hence, by the fundamental theorem of Young measures (see, e.g., [1]), 
for a suitable (not relabeled) subsequence, gk g in L X {T x [0,r]) where 

<j(M) = / k(A) - a(Du(:M))|di/( Xjt ) = 0 

JM( n + 1 ) Xm 

Here, we have used that by Lemma 11, we only have to integrate over and 

that there, by Lemma 12(i), <r(A) = cr(Du(x,t)). Since gk > 0, the convergence 
is strong. By the usual diagonal sequence argument, we can find a subsequence 
which converges strongly in iLCTxlO.oo]). □ 

This allows us now to finally prove the existence theorem for the heat flow 
of the energy E : 

Theorem 14. If the convex C 1 -function p satisfies the coercivity and growth con- 
ditions (14) and (15), then there exists a global weak solution of (16) for arbitrary 
initial data uq E W 1,p (T, S n ). 

Proof. It follows from Lemma 13, that 

a(Duk) — > cr(Du) in Lioc(T x [0, oo]) (44) 

By taking the wedge product of (18) with Uk, we find that 

0 = d t u k A Uk- div(p , (|£>u fc |) A u k ). (45) 

\Du k \ 

Passing to the limit in the weak form of (45), and using (33), (34) and (44) yields 

Du 

0 = d t u A u — div(p f (\Du\) ■ A u). (46) 

\Du\ 

As stated at the beginning of Section 2.2, \u\ — 1 almost everywhere on T x [0, oo[. 
Therefore, 

(d t u — div(a(Du)) — ucr(Du) : Du)mf = 0 (47) 

is an identity for arbitrary ^ E Cq°(T x [0, oo]). 

Let (j) E Cq°(T x [0,oo],M n+1 ). We test (46) by u/\<j) and use ^ = u-(f) in (47) 
(observe, that these are admissible test-functions) and subtract the resulting equa- 
tions from each another. Using the identity 

(j) = u(u •</>)— u A (u A (j)) 

this yields the weak form of (16). 



□ 




Heat Flow into Spheres for a Class of Energies 



63 



References 

[1] J.M. Ball: A version of the fundamental theorem for Young measures. In: Partial 
differential equations and continuum models of phase transitions: Proceedings of an 
NSF-CNRS joint seminar. Springer, 1989 

[2] F. Bethuel: The approximation problem for Sobolev maps between manifolds. Acta 
Math. 167 , 153-206 (1991) 

[3] F. Bethuel, X. Zheng: Density of smooth functions between two manifolds in Sobolev 
spaces. J. Funct. Anal. 80 , 60-75 (1988) 

[4] K.-C. Chang, W.-Y. Ding, R. Ye, R.: Finite-time blow-up of the heat flow of harmonic 
maps from surfaces. J. Differ. Geom. 36, No. 2, 507-515 (1992) 

[5] Y. Chen: The weak solutions to the evolution problem of harmonic maps. Math. Z. 
201 , 69-74 (1989) 

[6] Y. Chen, M.-C. Hong, N. Hungerbiihler: Heat flow of p-harmonic maps with values 
into spheres. Math. Z. 215 , 25-35 (1994) 

[7] P. Courilleau, F. Demengel: Heat flow for p-harmonic maps with values in the circle. 
Nonlinear Anal. 41, no. 5-6, Ser. A: Theory Methods, 689-700 (2000) 

[8] G. Dolzmann, N. Hungerbiihler, S. Muller: Non-linear elliptic systems with measure- 
valued right-hand side. Math. Z. 226 , 545-574 (1997) 

[9] F. Duzaar, M. Fuchs: On removable singularities of p-harmonic maps. Analyse non 
lineaire, 7 , No. 5 385-405 (1990) 

[10] F. Duzaar, M. Fuchs: Existence and regularity of functions which minimize certain 
energies in homotopy classes of mappings. Asymptotic Analysis 5, 129-144 (1991) 

[11] J. Eells, L. Lemaire: Another report on harmonic maps. Bull. London Math. Soc. 
20 , 385-524 (1988) 

[12] J. Eells, J.H. Sampson: Harmonic mappings of Riemannian manifolds. Am. J. Math. 
86, 109-169 (1964) 

[13] J. Eells, J.C. Wood: Restrictions on harmonic maps of surfaces. Topology 15, 263-266 
(1976) 

[14] A. Fardoun, R. Regbaoui: Equation de la chaleur pour les applications p-harmoniques 
entre varietes riemanniennes compactes. (French) C. R. Acad. Sci. Paris Ser. I Math. 
333, no. 11, 979-984 (2001) 

[15] A. Fardoun, R. Regbaoui: Heat flow for p-harmonic map with small initial data. To 
appear in Calc. Var. Partial Differential Equations 

[16] A. Fardoun, R. Regbaoui: Heat flow for p-harmonic maps between compact Rie- 
mannian manifolds. To appear in Indiana Univ. Math. J. 

[17] H. Federer: Geometric measure theory. Springer, 1969 

[18] M. Fuchs: Some regularity theorems for mappings which are stationary points of the 
p-energy functional. Analysis 9 , 127-143 (1989) 

[19] M. Fuchs: p-harmonic obstacle problems. Part III: Boundary regularity. Annali 
Mat. Pura Applicata 156 , 159-180 (1990) 

[20] R.S. Hamilton: Harmonic maps of manifolds with boundary. Lect. Notes Math. 471 , 
Springer, Berlin, 1975 




64 



N. Hungerbiihler 



[21] R. Hardt, F.H. Lin: Mappings minimizing the L p norm of the gradient. Comm. Pure 
and appl. Math 15 , 555-588 (1987) 

[22] N. Hungerbiihler: p-harmonic flow. Diss. Math. Wiss. ETH Zurich, Nr. 10740, 1994. 
Ref.: Michael Struwe; Korref.: Jurgen Moser, Zurich, 1994 

[23] N. Hungerbiihler: Global weak solutions of the p-harmonic flow into homogeneous 
spaces. Indiana Univ. Math. J. 45 / 1 , 275-288 (1996) 

[24] N. Hungerbiihler: Compactness properties of the p-harmonic flow into homogeneous 
spaces. Nonlinear Anal. 28/5, 793-798 (1997) 

[25] N. Hungerbiihler: Non- uniqueness for the p-harmonic flow. Canad. Math. Bull. 40/2, 
174-182 (1997) 

[26] N. Hungerbiihler: ra-harmonic flow. Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4) XXIV, 
4, 593-631 (1997) 

[27] N. Hungerbiihler: A refinement of Ball’s Theorem on Young measures. New York J. 
Math. 3 , 48-53 (1997) 

[28] N. Hungerbiihler: Young measures and nonlinear PDEs. Habilitationsschrift ETH 
Zurich, 2000 

[29] N. Hungerbiihler: Quasilinear parabolic systems in divergence form with weak mono- 
tonicity. Duke Math. J. 107 / 3 , 497-520 (2001) 

[30] L. Lemaire: Applications harmoniques de surfaces riemannienne. J. Diff. Geom. 13 , 
51-78 (1978) 

[31] X.G. Liu, S.H. Li: The p-harmonic heat flow with potential into homogeneous spaces. 
Acta Math. Sin. (Engl. Ser.) 18 , no. 1, 21-26 (2002) 

[32] X.G. Liu: A note on heat flow of p-harmonic mappings. (Chinese) Kexue Tongbao 
(Chinese) 42, no. 1, 15-18 (1997) 

[33] X.G. Liu: A remark on p-harmonic heat flows. Chinese Sci. Bull. 42, no. 6, 441-444 
(1997) 

[34] M. Misawa: On regularity for heat flows for p-harmonic maps. Proceedings of the 
Korea- Japan Partial Differential Equations Conference (Taejon, 1996), 14 pp., Lec- 
ture Notes Ser., 39 , Seoul Nat. Univ., Seoul, 1997 

[35] M. Misawa: Existence and regularity results for the gradient flow for p-harmonic 
maps. Electron. J. Differential Equations, No. 36, 17 pp. (1998) 

[36] M. Misawa: Existence and regularity results for the gradient flow for p-harmonic 
maps. Regularity, blowup and related properties of solutions to nonlinear evolution 
equations (Japanese). (Kyoto, 1997). Surikaisekikenkyusho Kokyuroku No. 1045 , 57- 
72 (1998) 

[37] M. Misawa: On the p-harmonic flow into spheres in the singular case. Nonlinear 
Anal. 50 , no. 4, Ser. A: Theory Methods, 485-494 (2002) 

[38] J. Sacks, K. Uhlenbeck: The existence of minimal immersions of 2-spheres. Ann. of 
Math. 113 , 1-24 (1981) 

[39] R.S. Schoen, K. Uhlenbeck: Approximation theorems for Sobolev mappings, preprint 

[40] P. Strzelecki: Regularity of p-harmonic maps from the p-dimensional ball into a 
sphere. Manuscr. Math. 82, No. 3-4, 407-415 (1994) 




Heat Flow into Spheres for a Class of Energies 



65 



[41] P. Strzelecki: Stationary p-harmonic maps into spheres. Janeczko, Stanislaw (ed.) et 
al., Singularities and differential equations. Proceedings of a symposium, Warsaw, 
Poland. Warsaw: Polish Academy of Sciences, Inst, of Mathematics, Banach Cent. 
Publ. 33, 383-393 (1996) 

[42] E. Zeidler: Nonlinear functional analysis and its applications, Vol. II/B: Nonlinear 
Monotone Operators. New York: Springer, 1990 

[43] C.Q. Zhou: Compactness properties of heat flows for weakly p-harmonic maps. (Chi- 
nese) Acta Math. Sinica 41 , no. 2, 327-336 (1998) 

[44] C.Q. Zhou, Y.D. Wang: Existence and nonuniqueness for the flow of p-harmonic 
maps. (Chinese) Acta Math. Sinica 41 , no. 3, 511-516 (1998) 

Norbert Hungerbiihler 

Department of Mathematics 

University of Fribourg 

P erodes 

CH-1700 Fribourg, Switzerland 




Progress in Nonlinear Differential Equations 
and Their Applications, Vol. 59, 67-80 
© 2004 Birkhauser Verlag Basel/Switzerland 



Singularity Models for the Ricci Flow: 
An Introductory Survey 

Dan Knopf 



Abstract. At the time this article is being written (May 2003) much of the 
mathematical world is waiting with intense interest to see the results of Gr- 
isha Perelman’s effort [22, 23] to resolve William Thurston’s Geometrization 
Conjecture [28] for closed 3-manifolds by completing the program [16] begun 
by Richard Hamilton. It is still too early to give an accurate and fair as- 
sessment of the full impact of Perelman’s work. But in order to aid the many 
mathematicians who may be inspired by that work to look more closely at the 
Ricci flow, this does seem like an appropriate time to write a brief and purely 
expository introduction to the topic, intended for the non-expert. Readers 
desiring more information are encouraged to read the more advanced survey 
articles [6] and [7], as well as to consult Hamilton’s and Perelman’s original 
papers. 



1. Heuristics 

There are two heuristic principles which are useful to keep in mind when one first 
studies the (unnormalized) Ricci flow, wherein one starts with a smooth Riemann- 
ian manifold (M n ,go) and evolves it by the equation 

r\ 

-g t 9 = - 2Rc, (1.1) 

where Rc denotes the Ricci tensor of the metric g. The first principle is that 
equation (1.1) is morally the heat equation for a Riemannian metric. The best 
way to see this is by writing the right-hand side of equation (1.1) in harmonic 
coordinates. Recall that a coordinate chart (p : U — > R n defined in a neighborhood 
U of a point Xo on a smooth Riemannian manifold (M. n ,g) is called harmonic if 
the coordinate functions {x fc :l</c<n}it induces are harmonic throughout U\ 

Ax k = 0. 

It follows from standard existence and regularity theory for elliptic PDE that a 
harmonic coordinate chart exists in a sufficiently small neighborhood U of any xo G 

Received by the editors Version 071103. 
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A4 n , and moreover that the metric enjoys optimal regularity in such coordinates. 
The feature of harmonic coordinates relevant to the Ricci flow is that the identity 

0 = -Ax k = -g ij ( - d l . - rt-i) X k = gVrtj 
\dx l dxJ 13 Ox 3 ) 13 

holds throughout U. Using this, it is not hard to see that 

- 2 Rij = - 2 RZj = 9 ke -^i (9a) + Qa (g-\ dg) , ( 1 . 2 ) 

where Q is quadratic in the inverse and first derivatives of g. In other words, the 
highest-order derivative term in —2 appears in a harmonic chart to be the 
Laplacian of the component of the metric, regarded as a scalar function in that 
chart. 

There are at least two reasons why this principle is heuristic and not entirely 
rigorous. One is that equation (1.2) is not tensorial; indeed, all covariant derivatives 
of the metric vanish identically. Another reason is that when one evolves a metric 
by (1.1), coordinates which are harmonic at time t cannot be expected to be so at 
time t + £, for any e > 0. 

Nonetheless, much of what this heuristic principle suggests is almost true. 
Equation (1.2) is a quasilinear parabolic equation. (It is quasilinear because the 
inverse of the unknown function g (t) multiplies the highest-order derivatives in the 
equation.) In fact, the Ricci flow is quasilinear and almost parabolic. Its failure 
to be strictly parabolic stems from the fact that the Ricci flow is defined entirely 
in terms of natural geometric quantities, and hence is invariant under the full 
diffeomorphism group. This invariance is a great advantage from the standpoint 
of geometry. Fortunately, it is only a minor inconvenience from the standpoint 
of pde, because Dennis DeTurck showed [10, 11] that the flow is equivalent to a 
strictly parabolic quasilinear equation, so that short-time existence and uniqueness 
follow from standard theory. 

The heuristic of thinking of the Ricci flow as a heat equation is useful in 
another important sense. The heat equation seeks to regularize its initial data, so 
one expects equation (1.1) to improve the metric, at least for a short time. In fact, 
initial bounds on the curvature of a metric induce subsequent a priori bounds on 
all derivatives of the curvature. These derivative bounds were derived in [3], [2], 
and [25, 26]; most known proofs use the maximum principle for parabolic equations 
in a familiar technique pioneered by Bernstein. 

On a more intuitive level, the heat equation heuristic leads one to expect that 
the Ricci flow will mimic the diffusion properties of the heat equation, and thus 
will try to make a metric more homogeneous and isotropic. As we will see below, 
this expectation too is often justified. 

The second heuristic principle is that we should expect the Ricci flow to 
develop singularities; in particular, it tells us that the first principle can be mis- 
leading if we take it too literally. One arrives at the second heuristic from the 
viewpoint of geometry when one observes that the normalized Ricci flow (defined 
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in (2.1) below) can converge only to an Einstein metric. Since most Riemannian 
manifolds of dimension n > 2 are not Einstein, one should expect something to 
go wrong. The second principle is also supported by the viewpoint of analysis. 
Indeed, equation (1.1) implies that the scalar curvature evolves by 

^R = Ai? + 2|Rc| 2 . (1.3) 

This is a reaction-diffusion equation: the reaction term 2 |Rc| 2 may be regarded 
as fighting against the diffusion term A R. By a standard estimate, equation (1.3) 
implies that 

^ R > A R+-R 2 . 
at n 

By the parabolic maximum principle, one can then compare solutions of (1.3) with 
solutions of the ode 

dp _ 2 ^ 2 
dt n ^ 

obtained by ignoring the diffusion term in (1.3). One concludes that 

o) -*(»-<«> 

for all t > to that the solution g (t) exists. Hence if the scalar curvature ever 
becomes everywhere positive, a finite-time singularity is inevitable. 

The second heuristic principle suggests correctly that a positive resolution of 
the Geometrization Conjecture depends on developing an adequate understanding 
of singularities. One must learn what the behaviour of a solution to the Ricci 
flow that becomes singular reveals about the topology of the underlying manifold. 
The reason why current research into geometrization concentrates on the analysis 
of singularities is because any solutions which remain nonsingular behave very 
nicely. Indeed, these are the solutions for which the diffusion term dominates, in 
accord with the first heuristic principle. Hence our informal survey begins with 
such solutions. 



2. Nonsingular solutions 



By rescaling space and time, one can modify the Ricci flow so that it preserves 
volume. This leads to the equation 



3 _ _ 2r 

g = — 2Rc +— g, 
ot n 



( 2 . 1 ) 



where 



r(t) = 



f M n Rdg 
Im 71 dp 



denotes the average scalar curvature. Equation (1.1) (the unnormalized Ricci flow) 
is in a sense the more natural pde, but equation (2.1) (the normalized Ricci flow) 
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can be more convenient for taking limits and for establishing certain convergence 
properties. 

For example, in proving Hamilton’s seminal result [14] that any manifold 
which admits a metric of positive Ricci curvature is in fact a space form, one first 
studies the solution of equation (1.1). One shows that the solution exists on a 
maximal time interval 0 < t < T < oo, that the curvature becomes unbounded as 
t /* T, and that the metric becomes nearly Einstein at points where the curvature 
is large. Then one derives a crucial estimate on the gradient of the scalar curva- 
ture. Together with the observation that the diameter of the solution is uniformly 
bounded, this allows one to show that the ratio R m ax /Rmin approaches 1 as the 
singularity time is approached. Then one rescales space and time, converting the 
solution g(t) of equation (1.1) into a solution g{t) of equation (2.1). In the final 
step, one shows that g (t) exists for all time and converges exponentially in every 
C k norm to a metric of constant sectional curvature. 

One way of understanding the role of the normalized flow in this proof is as 
follows: because the minimum curvature of the original solution approaches the 
maximum curvature as the maximum becomes large, the average curvature itself 
becomes large. Intuitively this means that when one rescales, the dilation term 
2r/n is large enough to keep the solution nonsingular. So in this case, the diffusion 
effect suggested by the first heuristic principle wins: it tames the singularity one 
expects from the second heuristic principle. On the other hand, if the curvature 
of a solution of (2.1) ever becomes large on a set of small volume, one would not 
expect the dilation term to be adequate to avoid a singularity. 

One says a solution (A i 3 ,g(t)) of the normalized flow (2.1) on a compact 
3-manifold is nonsingular if it exists for all positive time and satisfies a uniform 
curvature bound 

sup |Rm| < C < oo. 

A4 3 x [0,oo) 

In this case, A4 3 is geometrizable. This result was proved in [18]. The proof uses the 
Gromov- type convergence results in [17] as well as the advanced tensor maximum 
principle of [15], which lets one compare a tensor evolving by a PDE with a solution 
to a system of ODE. 

The results in [18] classify nonsingular solutions. The behaviours one ob- 
serves for such solutions turn out to be instructive when one later studies singular 
solutions. Hence it will be useful to recall some aspects of that classification here. 
If 

lim ( sup injfy,£) ) =0, (2.2) 

t_> °° \xeM 3 ) 

one says the solution exhibits collapse with bounded curvature. In this case, results 
of Cheeger- Gromov imply that M 3 admits an ^-structure, hence is a graph man- 
ifold, hence can be decomposed into a union of Seifert fiber-space pieces, all of 
which are known to be geometrizable. (An excellent survey is [24].) In case (2.2) 
does not hold, one can find £ > 0 and sequences of points Xj G A4 3 and times 
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tj — > oo such that 

inf(in j( ^))>e. (2.3) 

Then there exist diffeomorphisms tpj : M 3 — > M 3 such that the pointed sequence 
of solutions 

(M 3 ,gj ( t),Xj ) 

defined by 

9j (t) = (‘Pjff) (tj + t) 

converges locally smoothly in the pointed category to a limit solution 

(^oo) Qoo (t) 7%Oo) • 

If M^ is compact, it is necessarily diffeomorphic to M 3 . In this case, results 
in [18] prove that g ^ is a metric of constant sectional curvature K, hence that M 3 
is a space form. Moreover, one can in some cases prove a stronger convergence than 
was stated above. Indeed, if K > 0, then there exists tj large enough so that the 
Ricci curvature of g (tj) is strictly positive; it then follows from Hamilton’s original 
result [14] that the original solution converges exponentially fast to a space form. 
(In other words, one needed neither to modify by diffeomorphisms nor to restrict to 
a sequence of times.) If K = 0, the same statement holds by [13]. It is conjectured 
but not yet proven that a similar statement holds for the case K < 0. 

The remaining possibility is that M 3 ^ is noncompact, which is in a sense the 
most interesting case. A particular consequence of the results proved in [18] is that 
there exists a finite collection of finite-volume noncompact hyperbolic manifolds 
(7fy, hi) with the following properties. For each k E N, there exist truncations Hi,k 
of Hi along constant mean curvature tori of area less than 1/A:, a time tk, and 
embeddings (t) • H^k ^ Moo defined for all t > tk such that the images 
^Pi,k (t) (Hi,k) are mutually disjoint submanifolds for all t > tk and 

sup|K*fc(0 (9(t))~hi\\ ck{Hik) < 1/k 

t>t k 

in the norm of C k convergence on compact subsets of Hi : k- From the perspective 
of topology, the most important aspect of this case is that the fundamental group 
of each Hi injects into M 3 ^ under the maps c . This proves that M 3 ^ is Haken, 
hence is geometrizable by Thurston’s result [28]. In fact, one can show directly that 
the injectivity radius at any point in the complement Moo\^i^i,k (t) (Hi,k) of the 
hyperbolic pieces in the limit remains less than 1/k for all t > tk- In particular, 
this implies that the complement is a manifold-with-boundary which admits an 
^-structure, hence may be decomposed into a union of Seifert fiber-space pieces. 



3. Singular solutions 

Having seen that any manifold which admits a nonsingular solution of the Ricci 
flow is known to be geometrizable, one is led to ask the following question: What 
does a singularity of the flow reveal about the topology of the underlying manifold? 




72 



D. Knopf 



A powerful technique in geometric analysis for answering this question is to ‘blow- 
up’ the singularity, obtaining a limit flow whose properties should yield information 
about the geometry of the original manifold near the singularity just prior to its 
formation. 

To see how this is done, consider a solution (M 71 ,g (t)) of the unnormalized 
Ricci flow on a maximal time interval 0 < t < T. For simplicity, we shall discuss 
the case that T < oo, although an analogous theory is available for singularities 
which form in infinite time. In the case of a finite- time singularity, it follows from 
short-time existence results for the flow that 

lim ( sup |Rm(x,£)| I = oo. (3.1) 

ty'T J 

To develop intuition, consider what is perhaps the simplest example of a 
finite-time singularity. Let g can denote the canonical metric of constant sectional 
curvature K — 1 on the sphere S n . Consider a 1-parameter family g (t) of confor- 
mal metrics defined by 

g(t) = r(t) 2 g c an, (3.2) 

noting that g (t) has constant sectional curvature r(t)~ 2 . Observe that g(t) is a 
solution of the Ricci flow if and only if 

dv d 

— • <7can = ~ ^ Rc [g] = 2 Rc [^can] = 2 (n 1) • g c &n-> 

hence if and only if 

r ( t ) = y/r§-2(n-l)t = y/2 (n - l)y/T=t. (3.3) 

So if g (t) solves the Ricci flow, it must become singular at T = r\j (2 [n — 1)). 

In general, one says that a finite-time singularity of ( M n ,g(t )) is of Type I 
or fast-forming if it occurs at the natural rate suggested by this example, hence if 
there exists C < oo such that 

sup |Rm| • (T — t) < C < oo. 

M n x [0,T) 

On the other hand, one says the singularity is of Type Ila if 

sup |Rm| • (T — t) = oo. 

A4 n x [0,T) 

Solutions with this property are also called slowly- forming singularities. (This ter- 
minology may be somewhat nonintuitive initially, but is actually perfectly natural 
for reasons involving the sharpness of certain a priori estimates which guarantee 
short-time existence of the flow.) 

It is common to study a singularity by the technique of parabolic dilation. One 
chooses sequences of points Xj G M n , times tj E [0, T) increasing monotonically to 
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T, and dilating factors A j >0. Using these, one defines a sequence (Ad n , gj (t ) , Xj) 
of pointed solutions to the Ricci flow, where 

9j (t) = Xjg (tj + {) . 

Notice that each solution gj ( t ) exists on the time interval — Xjtj <t<Xj(T — tj). 

One chooses the sequences {xj}, {tj}, and {Ay } in order to get a good limit. 
A flat limit would not be desirable because it would not reveal enough about the 
geometry of the original solution. (Intuitively, getting a flat limit means one has 
dilated too much.) On the other hand, one wants uniform bounds on the curvatures 
of the metrics gj (0) in order to take advantage of the derivative estimates described 
in Section 1. (Bounds on the curvature of a solution to the Ricci flow imply bounds 
on all derivatives of the curvature.) Accordingly, one usually chooses Xj and Xj so 
that the curvature |Rm (xj, tj)\ is comparable to sup xe7V1 n |Rm (x, tj ) | and so that 

0 < c< XJ 1 |Rm(xj,£j)| < C < oo. (3.4) 

Recalling (3.1), one chooses the sequence of times so that sup^^n |Rm(x,{,-)| is 
comparable to the maximum curvature over a sufficiently large interval of earlier 
times. 

Having made suitable choices, one wants to prove that (A \ n ,gj (t) ,Xj) 
converges locally smoothly in the pointed category to a limit solution 
(•Moo?<7oo ( t ) ,Xoo) of the Ricci flow, called a singularity model (also called a ‘fi- 
nal time limit flow’ in the literature). As we shall see below, such solutions are 
typically very special. (For example, a singularity model formed from a Type Ha 
singularity will exist for all times — oo < £ < oo, which helps explain why such 
singularities are called ‘slowly forming’.) By developing an adequate understand- 
ing of these model solutions, one hopes to obtain local information about the 
geometry of the original sequence just prior to the formation of the singularity. In 
particular, the strategy for extracting topological information about a 3-manifold 
At 3 from a solution (A 4 3 ,g ( t )) of the Ricci flow that becomes singular is to per- 
form a geometric- topological surgery on Ad 3 just prior to the singularity in such 
a way that the maximum curvature of the solution is reduced by an amount large 
enough to permit the flow to be continued on the piece or pieces that remain af- 
ter the surgery. In the final step of this program, one needs to argue that only 
geometrically recognizable pieces will remain after finitely many surgeries. 

The main difficulty in obtaining convergence to a singularity model is ob- 
taining an adequate injectivity radius estimate for the sequence ( M n ,gj ( t ) , Xj), 
namely an estimate of the form 

inf 3 (inj (x,t)) > --= ° ■ 

xeM3 V su PxeM* l Rm {x, t) | 

Until recently, such estimates had to be derived by ad hoc means. For example, 
in Section 23 of [16], Hamilton proved an isoperimetric inequality that implies 
an injectivity radius estimate for appropriately chosen sequences of dilations ap- 
proaching a Type I singularity in dimension three. In Section 22 of the same paper, 




